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Abstract 

A pedagogical introduction to the heavy quark theory is given. It is explained 
that various expansions in the inverse heavy quark mass l/rng present a version of 
the Wilson operator product expansion in QCD. A systematic approach is developed 
and many practically interesting problems are considered. I show how the l/mq 
expansions can be built using the background field technique and how they work in 
particular applications. Interplay between perturbative and nonperturbative aspects 
of the heavy quark theory is discussed. 



*An extended version of the lectures given at Theoretical Advanced Study Institute QCD and 
Beyond, University of Colorado, Boulder, Colorado, June 1995. 



1 Lecture 1. Heavy Quark Symmetry 



The statement that Quantum Chromodynamics (QCD) is the theory of hadrons has 
become common place. It is a very strange theory, since many questions concern- 
ing dynamics of the quarks and gluons at large distances - however simple they 
might seem - remain unanswered or, at best, understood only at a qualitative level. 
Progress in the direction of the quantitative description of the hadronic properties 
is slow - every step bringing us closer to such a description is painfully difficult. 
At the same time new results, even modest, have a special weight for obvious rea- 
sons - QCD, unlike many other trendy theories in the modern high energy physics, 
definitely has a direct relation to Nature and will stay with us forever. 

Every hadron in a sense is built from quarks and/or gluons. I say "in a sense" 
because these are no ordinary building blocks. The number of degrees of freedom 
fluctuates and is not fixed; this we know for sure. At large distances we have to 
deal with a genuine strongly coupled field theory, and, as usual, the strong coupling 
creates complicated structures which can not be treated by perturbative methods. 
Then we feel helpless and are ready to use every opportunity, no matter where it 
comes from, if only it gives the slightest hope of getting a solid quantitative approach 
based on QCD. 

QCD has two faces, two components - hard and soft. The hard component is 
the realm of perturbative QCD. Not much will be said in these lectures about this 
aspect. Instead, we will concentrate on the soft component. Many years ago, at the 
dawn of the QCD era, it was noted that heavy quarks are, probably, the best 
probe of the soft component of the gluon fields out of all probes we have at our 
disposal. The developments we witnessed in recent years confirm this conclusion. 

The dynamics of soft degrees of freedom in QCD is the realm of non-perturbative 
phenomena. Having said this I hasten to add that there is an element of luck - 
transition from the perturbative regime to the non-perturbative one is very abrupt 
in QCD. In a sense the gauge coupling constant is abnormally small. I do not mean 
here the conventional logarithmic suppression of the running constant but, rather, 
the fact that 6, the first coefficient in the Gell-Mann-Low function, is numerically 
large. This fact allows us to forget, in the first approximation, about perturbative 
effects and focus on non-perturbative ones in a wide range of problems. It is more 
exact to say that we will concentrate on studying the soft degrees of freedom, but 
due to the fortunate circumstance of "abnormal" smallness of a<j(/x)/7r for as low 
normalization point as /i ~1 GeV, all effects due to the soft degrees of freedom are 
essentially non-perturbative. I will elucidate the precise meaning of this statement 
later. 

It would be great if we could just switch off - by adjusting some parameter - 
all hard processes in QCD without changing its soft component. Then we would be 
left with the confining dynamics in a clean and uncontaminated form; formulation 
of the theory would be much easier. The only parameter which might do the job is 
b. If we could tend h ^ oo with Aqcd fixed, the hard gluons would be suppressed 
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by powers of 1/6 while the soft component would presumably remain unaltered or 
almost unaltered. Unfortunately, nobody knows how to make the enhancement of b 
parametric. (The limit of the large number of colors, oo, does not work since, 

although b is definitely proportional to Nc in this case, the perturbative expansion 
for all planar graphs goes in Nc/b, not in 1/6 |^.) Therefore, we will have to rely on 
the numerical enhancement of b. In the first lectures I will merely assume that the 
hard gluon exchanges are non-existent. Later on, at the very end, we will return to 
this issue and will briefly discuss the impact of hard gluons. 

The purpose of these lectures is mainly pedagogical - the coverage of the topic is 
neither chronological nor comprehensive. Technically sophisticated issues and cal- 
culations are avoided whenever possible; instead I discuss particularly illuminating 
problems, in a simplified setting. The readers interested in specific advanced appli- 
cations (e.g. combining the l/mg expansions with the chiral perturbation theory 
1^) are referred to the original publications and the review papers summarizing 
a wealth of results obtained in the heavy quark theory after 1990. The presentation 
of the heavy quark theory below as a rule does not follow the standard pattern and 
is, rather, complementary with respect to the more traditional reviews We try 
to emphasize that the heavy quark theory and the heavy quark expansion is nothing 
else than a version of the Wilson operator product expansion (OPE) an aspect 
which usually remains fogged. 

1.1 Why heavy quarks? 

The quark-gluon dynamics is governed by the QCD Lagrangian 

^ = -\g%G^^, + ^ + ^ Q(^ jp - mQ)Q = 
^ q Q 

Aight + EO(^^-"^Q)Q (1-1) 
Q 

where is the gluon field strength tensor, the light quark fields (m, d and s) are 
generically denoted by q and are assumed, for simplicity, to be massless while the 
heavy quark fields are generically denoted by Q. To qualify as a heavy quark Q 
the corresponding mass term rriQ must be much larger than Aqcd- The charmed 
quark c can be called heavy only with some reservations and, in discussing the heavy 
quark theory, it is more appropriate to keep in mind b quarks. The hadrons to be 
considered are composed from one heavy quark Q, a light antiquark q, or diquark 
qq, and a gluon cloud which can also contain light quark-antiquark pairs. The role 
of the cloud is, of course, to keep all these objects together, in a colorless bound 
state which will be generically denoted by Hq. 

Quite naturally in the heavy quark theory, the gamma matrices used are those 
of the standard representation. 
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With these definitions of the gamma matrices the left-handed spinor has the form 

^/'L = (1 + 75)^- 

The hght component of Hq, its hght cloud, Q has a complicated structure - the 
soft modes of the light fields are strongly coupled and strongly fiuctuate. Basically, 
the only fact which we know for sure is that the light cloud is indeed light; typical 
frequencies are of order of Aqcd- One can try to visualize the light cloud as a soft 
medium. The heavy quark Q is then submerged in this medium. If the hard gluon 
exchanges are discarded the momentum which the heavy quark can borrow from 
the light cloud is of order of Aqcd, and the corresponding uncertainty in the energy 
of the heavy quark is of order Aq^jj^/mg. Since these quantities are much smaller 
than ttlq this means, in particular, that the heavy quark-antiquark pairs can not 
play a role. In other words, the field-theoretic (second-quantized) description of 
the heavy quark becomes redundant, and under the circumstances it is perfectly 
sufficient to treat one single heavy quark Q within quantum mechanics, which is 
infinitely simpler, of course, than any field theory. Moreover, one can systematically 
expand in l/rriQ. Thus, in the limit ttlq/Aqcd 00 the heavy quark component 
of Hq becomes easily manageable allowing one to use the heavy quark as a probe 
of the light cloud dynamics. The special advantages of this limit in QCD were first 
emphasized by Shuryak [^]. 

1.2 Descending Down 

In field theory one has to specify the normalization point fi where all operators are 
defined; in particular, the gauge coupling constant g and the quark mass ttiq are 
functions of fi. The original QCD Lagrangian ( |1 . 1| ) is formulated at very short dis- 
tances, or, which is the same, at a high normalization point /i = Mq where Mq is the 
mass of an ultraviolet regulator. In other words, the normalization point is assumed 
to be much higher than all mass scales in the theory, fi ^ mq. Constructing an 
effective theory intended for description of the low-energy properties of the heavy 
fiavor hadrons we must evolve the Lagrangian from the original high scale Mq down 
to a normalization point fi lying below the heavy quark masses tjiq. By evolving 
down I mean that we integrate out, step by step, all high-frequency modes in the 
theory thus calculating the Lagrangian describing dynamics of the soft modes, 
with characteristic frequencies less than fi. The hard (high-frequency) modes de- 
termine the coefficient functions in C{fi) while the contribution of the soft modes is 
hidden in the matrix elements of (an infinite set) of operators appearing in C{fi). 
This approach, which in the context of QCD was put forward by K. Wilson long ago, 
has become common. It is widely recognized and exploited in countless applications 
- from the ancient problem of the K meson decays to fresh trends in the lattice 

^In some papers devoted to the subject the hght cloud is referred to as 'brown muck'. I think 
it is absolutely unfair with respect to the soft components of the quark and gluon fields to call 
them 'brown muck' only because we are not smart enough to fully understand the corresponding 
dynamics. 
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calculations [0]. The peculiarity of the heavy quark theory is due to the fact that 
the in and out states we deal with contain heavy quarks. Therefore, although we do 
integrate out the field fluctuations with the frequencies down to /x the heavy quark 
fields themselves are not integrated out since we will be interested in physics in the 
sector with the Q charge 7^ 0. The effective Lagrangian C{ii) acts in this sector. 

If QCD was solved we could include in our explicit calculation of the effective 
Lagrangian all modes, descending down to = 0. The Lagrangian obtained in this 
way would be built in terms of the fields of physical mesons and baryons, not in 
terms of quarks and gluons, since the latter become irrelevant degrees of freedom in 
the infrared limit /i ^ 0. This Lagrangian would give us the full set of all conceivable 
amplitudes and would, thus, represent the final answer for the theory. There would 
be no need for any further calculations - one would just pick up the amplitude of 
interest and compare it with experimental data. 

This picture is quite Utopian, of course. The real QCD is not solved in the 
closed form, and in doing explicit calculations of the coefficients in the effective 
Lagrangian one can not put n = 0. The lower the value of fi the larger part of 
dynamics is accounted for in the explicit calculation. Therefore, we would like to 
have fi as low as possible; definitely /i <^ mq. The heavy quark can be treated as a 
non-relativistic object moving in the soft background field only provided the latter 
condition is met. On the other hand, to keep theoretical control over the explicit 
calculations of the coefficient functions we must stop at some /i ^ Aqcd, so that 
as(/i)/7r is still a sufficiently small expansion parameter. In practice this means that 
the best choice (which we will always stick to) is ~ several units times Aqcd- All 
coefficients in the effective Lagrangian obtained in this way will be functions of fi. 

Since fi is an auxiliary parameter predictions for physical quantities must be fi 
independent, of course. The fi dependence of the coefficients must be canceled by 
that coming from the physical matrix elements of the operators in C{fi). However, 
in calculating in the hard and soft domains (i.e. above fi and below fi) we make 
different approximations, so that the exact fi independence of the physical quantities 
can be lost. Since the transition from the hard to soft physics is very steep one may 
hope that our predictions will be very insensitive to the precise choice of fi provided 
that fi ~ several units times Aqcd- Below, if not stated to the contrary we will 
assume that the normalization point fi is chosen in this way. 

In descending from Mq down to fi the form of the Lagrangian (|1.1|) changes, and 
a series of operators of higher dimension appears. It is important that all these 
operators are Lorentz scalars. For instance, the heavy quark part of the Lagrangian 
takes the form 

Aeavy = E (^(^ ^ " "^q)Q + 7^Q{^ / 2)cr + ^ ^QTQqTql+O ( ^] 
Q [ ^^Q r,q^Q J \^Q/ 

(1.3) 

where cg and dq^ are coefficient functions, Gfj,^, = gG'^^J,"' and is the color gen- 
erator, (Trt"t^ = 5°"^ /2)] below we will often use the short-hand notation iaG = 
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ia^jjG^i, = i'jfj.'juGniy. The sum over the hght quark flavors is shown exphcitly 
as well as the sum over possible structures F of the four-fermion operators. All 
masses and couphngs, as well as the coefficient functions cq and d^^\ depend on 
the normalization point. For example, the coefficient cg in the leading logarithmic 
approximation can be written as 

where Uf is the number of the light flavors. The power —3/6 was flrst calculated in 
Ref. 0. In Sect. 5.3 I will explain how to derive Eq. (|1.4|). 

The operators of dimension flve and higher in Eq. ( p..3|) are due to the contribu- 
tion of hard gluons, with offshellness from /i up to Mq. Since we agreed that in this 
lecture we will ignore the existence of such gluons, we will forget about these oper- 
ators for the time being. Does this mean that what remains from the Lagrangian 
( p..3|) contains no l/mq terms? 

The answer to this question is negative. The l/mq expansion is generated by 
the flrst ("tree-level") term in the Lagrangian (p..3|). 



Cl,., = Qif-mQ)Q. (1.5) 

Although the fleld Q in this Lagrangian is normalized at a low point fi the fleld Q 
carries a hidden large parameter, mq; isolating this parameter opens the way to the 
l/mq expansion. Indeed, the interaction of the heavy quark with the light degrees 
of freedom enters through = iD^, where 

The background gluon fleld is weak if measured in the scale mg, which means, of 
course, that there is a large "mechanical" part in the x dependence of Q(x), known 
from the very beginning 0, 

Q[x) = e-^"^'3*g(x) (1.6) 

where Q{x) is a "rescaled" bispinor fleld which, in the leading approximation, carries 
no information about the heavy quark mass. It describes a residual motion of the 
heavy quark inside the heavy hadron [|lOl with typical momenta of order Aqcd- 
Remnants of the heavy quark mass appear in Q only at the level of l/mg corrections. 
Equation (L6) is written in the rest frame of Hq. In the arbitrary frame one 



singles out the factor exp{—imqv^x^) where the four-velocity of the heavy hadron. 

The covariant momentum operator acting on the original flled Q, when acting 
on the rescaled fleld Q, is substituted by the operator mqv^ + vr^, 

zD^Qix) = e-'^Q^-^- {mqv^ + iD^) Q{x) = e"*'"'?''"^- {mqv^ + vr^) Q{x) . (1.7) 



5 



Below we will consistently use different letters, and vr^ for the momentum op- 
erators iD^ acting on Q and Q, respectively. If not stated to the contrary, we will 
use the rescaled field Q, omitting the tilde in all expressions where there is no risk 
of confusion ^ In the local colorless operators bilinear in the heavy quark fields it 
does not matter whether the original field Q or the rescaled one is used, since, say, 

qq = Qq, qv^,q = Qtt^q,... 

and so on. Using these distinct notations for the momentum operator is convenient 
since all expressions written in terms of vr^ and Q do not contain implicitly the large 
parameter mg. 

I pause here to make a reservation. The rescaled field Q is a four-component'Diia.c 
bispinor, not a two component non-relativistic spinor which is usually introduced 
in the heavy quark effective theory (HQET) DTof. HQET is a formalism invented in 
the very beginning of the 90 's |Ty] which is very often used in connection with the 
heavy quark physics @]. It is convenient in a range of problems but can be quite 
misleading in some other problems. I prefer to discuss the heavy quark expansions 
directly and systematically in full QCD in the framework of the Wilson OPE. In 
many instances the careful reader will certainly recognize a significant overlap, but 
the Wilson language, being more general, seems to give a better understanding and 
command over the l/mg expansions. Moreover, some issues can not be addressed 
in the framework of HQET at all. 

The Dirac equation (J/ — mQ)Q = in terms of the rescaled field can be written 
as follows: 

1 - 7o ^ i 



and 



T^oQ = ^ Q- 1-9 



The last equation is actually the squared Dirac equation. 



(f + mQ){f-mQ)Q = -^{V'+'-aG-ml)Q = 0. 



2mQ 2mg 

In deriving Eq. (|1.9|) we used the fact that 

[V„V,] = [n„n^]=tgG';,J\ (1.10) 

Armed with this knowledge one can easily obtain the l/mg expansion of C^eavyy 
up to terms l/mg. 



■^Whenever one sees an expression containing tt's one may be sure that it refers to the rescaled 
fields Q even if the tildes are not written out explicitly 
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i,(-,5i)H-2..ix,)](..||)i±-«.o(4). (i.n, 



where a denote the Pauh matrices and 



(jcaf = 7i'^ + aB . 



E and B denote the background chromoelectric and chromomagnetic fields, respec- 
tively. The coupling constant g and the color matrix are included in the definition 
of these fields. The derivation of this Lagrangian is a good home exercise. I encour- 
age everyone to obtain Eq. ( 1 . 1 1| ) by using the commutation relation ( |1.10 ) and 



the properties of the gamma matrices. Those who will have problems with getting 
Eq. ( p.. Ill) should consult Chapter 4 of Bjorken and Drell |jll| or Sect. 33 of the 



Landau-Lifshitz course []T2| from where this Lagrangian follows immediately. It is 
worth noting that 

Cl,.y^^^ino-nQ)cp (1.12) 

where 

and Tig is a non-relativistic Hamiltonian, through second order in I/itlq, 



well-known (in the Abelian case) from the text-book expressions |TT|, |T2[. Equation 



( p..l3| ) is merely the Foldy-Wouthuysen transformation which is necessary to keep 



the term linear in ttq in its canonic form. 

1.3 niQ oo; The heavy quark symmetry 

Let us first neglect all l/rriQ corrections altogether. In this limit tuq drops out from 
ro 

■^heavy ' 



,1 + 70 
2 

This expression takes place in the rest frame of Hq; in the arbitrary frame 



^Lvy = 0^vrog. (1.15) 



Cl^^y = Q^TT.v.Q . (1.16) 

In the limit tuq oo the masses of all Q-containing hadrons become equal to that 
of the heavy quark Q, 

Mhq =mQ + O(Aqcd) • 

The mass splittings between different hadrons are generically of order Aqcd ^ 
rriQ. Soon, we will relate these mass splittings to the expectation values of certain 
operators. 
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The assertion that all Q-containing hadrons are degenerate to the zeroth order 
in niQ is trivial. This "degeneracy" by no means implies that the internal structure 
of all Q-containing hadrons is the same. A little less trivial is the fact that there 
exist hadrons whose masses are degenerate to much better accuracy, 0{mQ^), and 
whose internal structure is, indeed, identical in the limit mq —>■ oo. 

Since all effects due to the heavy quark spin are, obviously, proportional to 
l/rriQ, in this limit the heavy quark spin becomes irrelevant, see Eqs. ( |1 . 1 1|) , ( |1.16| ). 
Correspondingly, there emerges a symmetry between the states which differ only by 
the spin orientation of the heavy quark. The pseudoscalar and vector mesons of 
the type B and B* (both are the ground state S wave mesons) present an example 
of such spin family. In the limit tuq oo their masses must be degenerate up to 
terms O^ttiq^), and the light clouds of B and B* coincide. If there is more than one 
heavy quark, say Qi and Q2, the theory is symmetric with respect to the interchange 
Qi ^ Q2 even if their masses are not close to each other (in physical applications 
we, of course, keep in mind b and c). Indeed, the heavy quark Qi plays the role of 
the static force center inside Hq.; the light cloud is flavor-blind and does not notice 
the substitution of Qi by Q2 provided that the four-velocities of both quarks are the 
same. Notice that at this level the four-velocity of the heavy quark coincides with 
that of the heavy hadron. (Only when higher order corrections in l/mq are taken 
into account the difference between the four-velocities becomes important and the 
symmetry Qi ^ Q2 is violated. At the level of l/mq also the spin symmetry is 
not valid any more.) If the hard gluon effects are neglected the interaction with the 
light cloud can not change the heavy quark four-velocity; therefore, this quantity is 
conserved in the strong interactions [1^. (This conservation is, of course, destroyed 
by the hard gluons which can easily carry away a finite fraction of the heavy quark 
momentum.) 

The symmetry connecting Qi and Q2 emerges in the limit mg^ ^ —>■ 00 even if 
the masses of the heavy quarks are not close to each other. What is important is 
that both must be much larger than Aqcd- We encounter here a situation which is 
conceptually close to the problem of the isotopic symmetry of the strong interactions. 
Everybody knows that the strong amplitudes are isotopically invariant with the 
accuracy up to a few percent, and, at the same time, the masses of the d and u 
quarks are not too close to each other, md/mu ~ 2. It is not the proximity of these 
masses which counts, but the fact that the both masses are much less than the QCD 
scale Aqcd- 

Usually the existence of an internal symmetry implies a degeneracy of the spec- 
trum. For instance, the isotopic symmetry mentioned above, apart from certain 
relations between the scattering amplitudes, predicts that the proton and neutron 
masses are the same, up to small corrections due to the symmetry breaking effects. 
The heavy quark symmetry does not manifest itself as a degeneracy in the spectrum 
- the D and B masses are very far from each other. One has to subtract the me- 
chanical part of the heavy quark mass in order to see that all dynamical parameters 
are insensitive to the substitution Qi Q2 in the limit mq^ ^ 00 |]TB[. Perhaps, 
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this is the reason why it was discovered so late. 

To elucidate the issue of the heavy quark symmetry let us consider a practical 
problem, semileptonic decay of the B meson induced by the weak b ^ c transition. 
The initial B meson decays into an electron-neutrino pair plus the D meson. Since 
we do not now discuss the l/rriQ corrections we may make no distinction between 
the four- velocities of the quark Q and the hadron Hq, and between their masses. 
Assume that the B meson is at rest. Furthermore, let us assume that the four- 
momentum q carried away by the lepton pair is maximal, = {Mb — M^y. This 
means that the D meson produced is also at rest - the hadronic system experiences 
no recoil. The corresponding regime is sometimes called the point of zero recoil. 

In this regime the B ^ D transition form factor is exactly unity! More exactly, 

{D\c-fob\B) = (2Mb2Mb)^/2 x unity (at zero recoil) (1.17) 

where the square root factors are due to the relativistic normalization of our ampli- 
tudes. By the same token 

{D*\cja5b\B) = z(2Mb2Md)^/^A* x unity (at zero recoil) (1.18) 

where D* is the polarization vector of D*. As well-known, the exact relations of this 
type always reflect an underlying symmetry. They can never emerge accidentally 
because only a sjTiimetry can protect the form factors from renormalizations. 

It is very easy to understand why Eqs. ( |1.17| ) and ( |1.18D take place. Indeed, 



the space-time picture is very transparent. The b quark at rest is surrounded by its 
light cloud, the latter being the eigenstate of the problem of color interaction with 
a static force center. At time zero the weak current instantaneously substitutes the 
b quark by c; the charmed quark is also at rest, and since the color interactions 
are flavor-blind the same light cloud continues to be the eigenstate, this time with 
the c quark as the static center. If instead of the field-theoretic light cloud we had 
a quantum-mechanical problem one could say that the overlap integral for these 
identical wave functions is 1. The light cloud will feel the substitution b ^ c only 
to the extent the heavy quark momentum inside the heavy meson does not vanish 
exactly - this effect is, of course, suppressed by powers of l/mq. As we will see later 
corrections in the right-hand side of Eqs. ( |1 . 1 71) and ( |1.18D are actually of order 



l/mg; there are no linear corrections in I/uiq. In the B —>■ D* transition generated 
by the axial- vector current the current, additionally, changes the orientation of the 
heavy quark spin. As was already mentioned, all effects related to the heavy quark 
spin are suppressed by l/mg; D and D* are in the same multiplet, and the B D* 
transition is governed by the same symmetry. This symmetry allows one to rotate 
arbitrarily four states, 

b spin up, b spin down, c spin up, c spin down; 

therefore, we obviously deal here with an SU(4) invariance. 
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The symmetry relations ( |1 . 1 T\ ) and ( |1.18| ) were first derived in Refs. |T5 . 



Shortly after it was realized [jT6| that the actual symmetry is much stronger - the 
SU(4) invariance takes place for any given value of v^, the four-velocity of the re- 
coiling c quark, not necessarily at the point of zero recoil or close to it. Thus, many 
different form factors connecting [B, B*) and {D, D*) can be expressed in terms 
of one function depending only on the velocity of the recoiling hadron (in the rest 
frame of the decaying hadron). The universal form factor is called the Isgur-Wise 
function. 

1.4 The Isgur-Wise function 

Now we are finally ready to discuss a very elegant observation due to Isgur and Wise 



T6| . Let us consider now the amplitudes induced by the transition cTh off the zero 
recoil point. Here F is any Lorentz matrix; of special interest are, of course, the 
vector and the axial-vector cases, 

r = 7^ or 7^75; 

the weak decays of the B meson are induced by the V ~ A currents. The physically 
measurable amplitudes are {D\cTh\B) and {D*\(Th\B)] for completeness one can 
also consider the amplitudes of the type {D\(Th\B*) or {D\cTh\B*) , or {B\bTh\B) - 
this adds nothing new. The four-velocity of the particle Hq is defined as 

y = ^Ilsk. (1.19) 

the four-velocities of the initial particles will be denoted by v while those of the 
final particles by v'. It is obvious that = 1, and, additionally, in the rest frame 
V = {1, 0, 0, 0}. In the most general case the amplitude (D|c7^6|_B) can be expressed 
in terms of two form factors, the amplitude (D*|c7^756|-B) in terms of three form 
factors and the amplitude (D*|c7^6|i?) in terms of one form factor. The heavy 
quark symmetry tells us that in the limit ttiq — > 00 these six functions, a priori 
independent, reduce to one and the same function which depends only on the scalar 
product vv'. Specifically, 



{D\c^,b\B) = JMbMo {v + v')^ ay) ; (1-20) 



{D*\c^,j,b\B) = JMbMdi DlH+vv') - (/^X)< i{y) ; (1.21 



{D*\ci,h\B) = \JMbMd [-e^.xaDlv'^v^\ ^y) . (1.22) 

where 

y = vv' 

and C,{y) is the Isgur-Wise function, is the polarization vector of D*. The Isgur- 



Wise function is independent of the heavy quark masses. The square root ^MbMd 
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reflects the relativistic normalization of the states. The symmetry relations ( |1.17| ) 
and (|1.18| ) imply that the normalization of the Isgur-Wise function at zero recoil is 
fixed, 

e(y=l) = l. (1.23) 



Perhaps, it is worth noting that the phases in Eqs. ( |1.21| ) and (|1.22|) differ from 



what you might see in the literature. They are, of course, a matter of convention 
and reflect the definition of the states. The definition I follow is in accord with the 
standard relativistic convention, see Eq. (|1.25| ) and Sect. 3.5. 

The fact that a large set of form factors degenerate into a single function de- 
pending only on y might seem a miracle; but after the assertion is made, with the 
knowledge you already have, it should be not difficult to understand why it happens. 
Indeed, let us turn again to the space-time picture described above. A h quark at 
rest, surrounded by the light cloud, instantaneously converts into a c quark. This 
time the four-momentum carried away by the lepton pair is not maximal; therefore, 
the c quark is not at rest. This force center flies away with the velocity v' . But the 
light cloud stays intact. So, the question is: "what is the amplitude for the flying c 
quark and the cloud at rest to form a D or D* meson?" We can look at this process 
in another way. After the b ^ c transition happened let us proceed to the rest frame 
of c. In this reference frame the c quark produced is at rest, but the cloud, as a 
whole, moves away with the velocity —v'. It is clear that this system - the static 
charmed force center plus a moving light cloud - has a projection on D or D* . The 
amplitude per se, with the kinematic structures excluded, can depend only on {ifl - 
there is no preferred orientation in the space, and the direction of is irrelevant. 
Using covariant notations one can say that the amplitude depends only on vv' since 



in the B rest frame vv' = y (1 + v''^). There is simply no place for the dependence 
on the heavy quark masses, apart from the overall normalization factors appearing 
because we stick to the relativistic normalization of states ^ 

Since the heavy quark spin is irrelevant in the limit ttiq — *^ cxo to warm up let 
us consider a toy model where the heavy quarks are deprived of their spins from 
the very beginning. In other words, I replace the genuine spin-1/2 heavy quarks of 
QCD by spin-0 color triplets with the same mass. We will turn to this toy model 
more than once below. 

In QCD, B and B* form a multiplet which includes 4 states: the total angular 
momentum of the light cloud (1/2) combines with the heavy quark spin (1/2) to 
produce either spin-0 state (B) or three spin-1 states [B*). In our toy model the 
analog of this ground-state multiplet is obviously a baryon of spin 1/2; let us denote 
the corresponding field by Nq, where a is the spinorial index. The current generating 
the transition A*";, —>■ Nc has the form c^b, where the fields b and c are assumed to 
be scalar now. At first sight the amplitude {Nc\c^b\Niy) might contain four different 



^Warning: an additional dependence on the heavy quark masses may emerge if we include the 
hard gluon exchanges neglected so far. For details see Ref. 
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kinematic structures, 



where I list only P-even structures, of course. On mass shell they all reduce to the 
first one, however; for instance, i/Nb = Nb- Hence 

{mc^b\Nb) = NcNb^y) . 

Returning to real QCD what remains to be done is to work out the consequences 
of spin. The most concise general formula can be written in terms of the matrices 

M = 5(z75) + 5^7/. and M' = D*{i^,) + D;^^ , (1.24) 

where 5^ and D^^ are the polarization vectors of B* and D*, respectively. The 
(heavy quark) spin independence of the strong interactions at itlq oo manifests 
itself in the fact that the couplings of the ground state pseudoscalar to i'j^ and the 
ground state vector to 7^ are the same, see Sect. 3.4. Now, the whole set of the 
transition amplitudes can be expressed by one compact cute formula Q , 

^ -.{H,{v')\crb\Hb{v)) = l-TT{M'^^r^M}ay = vv'). (1.25) 



Completing the trace we recover Eqs. ( |1.2q ) - ( |1.22| ). 



Equation (|L2|) can be derived in many different ways. Originally it was obtained 
in Ref. (see also 10). In Sect. 3.5 we will discuss one of the possible derivations 
- perhaps, not the simplest, but very instructive. Before we will be able to do that 
it is necessary to make a digression and study some elements of the background field 
technique. 

Equations ( |1.20[ ) - ( |1.22D are valid not only in the space-like domain (the form 



factor kinematics) but also in the time-like domain. The latter assertion calls for an 
immediate reservation, though. The heavy quark symmetry implies that Mb = Mb* ■ 
In the real world this equality is not exact: the heavy quark symmetry is violated 
by small l/mq terms. This small violation can be strongly enhanced in the near 
threshold domain, E ^ 2Mb ,where the symmetry breaking parameter turns out to 
be of order one [|19|. Indeed, let us consider a kinematical point above the threshold 



of the BB production but below BB*. In this domain all form factors describing 
three amplitudes 

{BB\J^\0), {BBV,\0), (5*5*1 J^|0), 

(here is some heavy quark current, say, = by^b) have imaginary parts associ- 
ated with the normal thresholds due to the intermediate state BB. On the other 
hand, there is no contribution to the imaginary part from the intermediate state BB* 



^Strictly speaking, for the outgoing particles one must use A4 — 7o(A^')Wo- With our conven- 
tions, however, Ad' = A4'. 
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and B*B* . In the pseudoscalar meson B the spin of the heavy quark Q is rigidly 
correlated with that of the light cloud. Hence, the spin independence of the heavy 
quark interaction is totally lost in the imaginary part in this point. In particular, 
in the amplitude {B* B*\J a kinematic structure forbidden by the Isgur-Wise 
formula appears. An even more pronounced effect of the heavy quark symmetry 
violation takes place in the anomalous thresholds generated by the pion exchange 
which can start parametrically much below the the normal thresholds depending on 



the interplay between M^* — and the pion mass ||2^ 



1.5 The mass formula 

To complete our first encounter with the basics of the heavy quark theory we will 
now derive a l/mq expansion for the masses of the Q-containing hadrons. 

It is intuitively clear that the heavy hadron mass can be expanded in terms of 
that of the heavy quark as follows 

MHQ=mQ + A + 0{mQ') (1.26) 

where A is a constant, of order of Aqcd, which depends on the light quark content 
and the quantum numbers of Hq but is independent of ttlq. (It first appeared in Ref. 
2T[. Later on we will see that this expression is not as trivial as it might naively 



seem and requires thoughtful definitions of all parameters involved. In particular, 
since the quarks are never observed as isolated objects, one may ask what the quark 
mass TTiQ actually means. In due time we will return to this question, of course. For 
the time being we agreed to disregard hard gluon exchanges; then mq is just the 
mass parameter in the Lagrangian ( |1.1| ). 

Formally Eq. (|1.26|) can be most easily derived by analyzing the trace of the 



energy-momentum tensor in QCD, 



9,, = mgQQ + ^G%G% (1.27) 



where l3{as) is the Gell-Mann-Low function. For simplicity we assume the light 
quarks to be massless; introduction of the light quark masses changes only technical 
details at intermediate stages of our analysis. If the mass term of the light quarks 
is set equal to zero the light quark fields do not appear explicitly in the trace of the 
energy-momentum tensor. The expression ( |1.27| ) contains two terms: the first one 
is a mechanical part while the second term is the famous trace anomaly of QCD ||2^ 
(for a review see e.g. Ref. [P^] ). 

Furthermore, as well-known, for any given one-particle state the expectation 
value of the trace of the energy-momentum tensor reduces to the mass of the state. 
Then, the hadron mass can be expressed in terms of two expectation values, 

- ^J^'^^'-'^^QIH,) + ^JH,l^-^G^lH,) (1.28) 
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where the relativistic normahzation of the states is imphed 



{Hq\Hq) = 2Mn^V 

in the rest frame; V is the normahzation volume. We wiU always use only the rel- 
ativistic normalization of states which will routinely result in the factors {2Mhq)~^ 
in all expressions. 

Let us discuss the expectation values of the operators in Eq. (|1.28|) in turn. 



The first one is explicitly proportional to mq. To be more quantitative we must 
determine the matrix element of the heavy quark density QQ. To this end it is 
convenient to use an argument suggested in Ref. [^| which will show us that the 



expectation value of QQ is very close to unity; as a matter of fact, with our present 
accuracy it is just equal to unity. The second expectation value reduces to A. 

Indeed, in the rest frame of Hq a typical momentum of Q is of order Aqcd, i-e. 
the heavy quark is very slow. This means that the lower components of the bispinor 
field Q are small compared to the upper ones and, hence, the scalar density of the 
heavy quark is close to its vector charge, QQ ~ QloQ- The difference is only due 
to the lower components. The vector charge, however, just measures the number of 
the heavy quarks inside Hq\ therefore, its matrix element is exactly unity. 

It is instructive to do the simple derivation outlined above in some detail. Com- 
bining the equations of motion, ( |1.8| ) and ( |1.9D , it is easy to get that 

which implies, in turn, 

QQ = QloQ - ttVO (tt^ + (tB)Q + higher orders (1.30) 



where B is the chromomagnetic field, Bi = eijkGjk- Equation (|1.30|) is the desired 
result demonstrating that 

,^{Hq\QQ\Hq) = -^{Hq\Q^oQ\Hq) + OiniQ') = 1 + OiniQ') . (1.31) 

The matrix element of the vector charge (appropriately normalized) is set equal to 
unity, as was discussed above. 

This digression has been undertaken merely to familiarize the reader with the 
basics of the l/mq expansion in QCD. As our understanding progresses the level 
of the explanatory remarks will be reduced so that in the subsequent lectures many 
derivations of a more technical nature will be suggested as an exercise. 

Thus, we have established that the first expectation value in Eq. (|1.28| ) produces 
TTiQ in the expansion for the heavy hadron mass. The second expectation value 
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which also has the dimension of mass obviously does not scale with mg in the limit 
TTiQ — i> cx), so one can define ^ 

A = ^(Hq\^G'\HqU,^^^. (1.32) 

Thus, the parameter A of the heavy quark theory is, in a sense, similar to the gluon 
condensate ^T\. The latter is the expectation value of the same gluon operator over 



the vacuum state. In the case of A the gluon operator is averaged over the lowest 
state of the system with the given (unit) value of the heavy quark charge. The lowest 
state is, of course, the ground state pseudoscalar meson, B. Generally speaking, Hq 
can be any Q-containing hadron. B mesons are most interesting from the point 
of view of applications; of practical interest also are Q-containing baryons which 
are the lowest-lying states in the given channel with the baryon quantum numbers. 
Therefore, strictly speaking, unlike the gluon condensate, there exist many different 
parameters A, one for every channel considered. Usually we will tacitly assume that 
A is defined with respect to the B mesons. 
Both expectation values. 



{HQ\mQQQ\HQ) and -^{Hq\^G'\Hq) 



have l/rriQ corrections which show up at the level 0{mQ^) in Eq. ( |1.26| ). Later on 
we will derive the expansion for Mhq which takes into account these terms 0{mQ^). 

The l/mq corrections in the expectation value of the gluon anomaly are due to 
the fact that in our approach the states \Hq) are physical heavy fiavor states, rather 
than the asymptotic states corresponding to mq = oo which are usually considered 
within HQET. Instead of working with these fictitious states I prefer to explicitly 
keep track of all 1 / mq corrections, both in the operators and in the definition of the 
states, appealing directly to the Wilsonean operator product expansion. 



^This expression relating A to the expectation value of the gluon anomaly operator was obtained 
in Ref. P5|. Some subtleties left aside in the derivation presented here are discussed in detail in 



this paper. It is instructive to compare Eq. ( 1.32| ) with a similar expression for the nucleon mass, 



known from ancient times 
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2 Lecture 2. Basics of the Background Field Tech- 
nique 

The essence of our approach is separation of all momenta into two classes - hard 
and soft. For the time being we will continue to pretend that the role of the gluon 
degrees of freedom reduces to a soft gluon medium. This is an ideal situation where 
the gluons can be treated as a background field. A powerful method allowing one 
to put calculations in the background fields on an industrial basis was developed 
by Schwinger in electrodynamics many years ago. In the eighties it was adapted 
to QCD. We will be unable to submerge in all details of this technique, and will, 



rather, present some basic elements in particular examples. The review paper ||28 
is recommended for further education. This lecture will be rather technical - its 
primary goal is to teach how the heavy quark mass expansions can be constructed 
in a systematic way in different problems. 

The starting point of the method is decomposition of fields into two parts - the 
quantum part and the background one. The propagation of quanta is described by 
the correlation functions of the quantum part of the fields considered in the external 
field. Later on the external field is to be considered as a fluctuating field of the light 
cloud, but this stage need not concern us at the moment. 

Let us start with a brief review of the Schwinger method, as it can be applied in 
QCD. We introduce the coordinate and momentum operators, and p^, respec- 
tively, = ig^u, [p^,Pu] = = 0. Moreover, introduce a formal set of 
states \x) which are the eigenstates of the coordinate operator X^, 

Xf,\x) =x^\x). (2.1) 

Please, note that \x) has nothing to do with the field-theoretic eigenstates, e.g. \Hq). 
To emphasize this fact the use the regular bracket ) in the notation instead of the 
angle one, which is reserved for the field-theoretic eigenstates. 

Then define the covariant momentum operator satisfying the following com- 
mutation relations 

[V,, X,] = tg^,, [V^, V,] = igf^Gl,, (2.2) 

where t"" are the generators of the color group, G^^, is the external field. 

The algebra ( p. 21) is the basic tool of the Schwinger formalism. We will expand 
the Green functions in the background field, and in each order of the expansion we 
will need to use only this algebra. 

In the coordinate basis acts as a covariant derivative, namely 

{y\V,\x) = (^z^ + g^A^^ix)^ 6{x - y) (2.3) 

if 

{y\x)=6{x-y). (2.4) 
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Now we can write formal expressions for the Green functions. For instance, for 
the quark Green function (mass nig) describing propagation from the point to the 
point X we have 

Six,0) = ix\-^\0) . (2.5) 
Y -rUg 

Eq. ( |2.5| ), rather obvious by itself, is readily verified by applying the Dirac 
operator to both sides of Eq. (|2.5| ). Furthermore, it can be identically rewritten as 
follows 

g(^.0)^M(y + ^.) p,_,„,/(^,^^)^^^^^J O), (2.6) 

where 

Please, note that the ordering is important here since does not commute with 

Vy of G^y. 

If we are aimed at calculating the coefficient functions in the Born approximation 
we need nothing else - Eq. ( |2.(j| ) is just systematically expanded in powers of the 
background field by using the commutation relations ( p.2| ). 

Observe that one can always shift by a c-number vector due to the fact that 

^igX^^^-^,X _ ^ (2.7) 

Hence, the Fourier transformed propagator reduces to 

dW«^5(x,0) = / d^xix\e''i^—^ e-*'?^|0) = / d'^xix\- \ 10). 

J If -niq J lf+ q-niq 

This simple trick allows one to readily develop the expansion sought for. Indeed, 
assume that q is large (hard momentum) and V represents soft modes and is small 
in this sense. Then we can expand in P, 

rfVe'''"^(x, 0) ^ — ^ -^f-^ + ... (2.8) 

q— nig q— nig q— nig 

Next, we transpose V to the right-most (left-most) position and act on the states 
using the equations of motion. 

It may seem that so far we got almost nothing compared to the standard Feyn- 
man graph calculations. Let us demonstrate the efficiency of the background field 
technique in a few examples. 



2.1 Inclusive decay of the heavy quark — toy model 

One of the most important practical problems in the heavy quark theory is the 
description of the inclusive decays of heavy fiavors. The semileptonic and radiative 
decays of the B mesons B X^v and B X^'y are particular examples. Both 
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are two key elements of the ongoing experimental efforts, in quest of new physics. 
Needless to say that a reliable QCD-based theory of such decays is badly needed. 
In this section we start discussing basics of such a theory. 

Since this is our first exercise, for pedagogical reasons, it seems reasonable to 
"peel off" all inessential technicalities, like the quark spins, and resort to a simplified 
model. In this toy model we will consider the inclusive decay of a spinless heavy 
quark into a spinless lighter quark plus a photon. Of course, our photon is also a toy 
photon. We will assume it to be scalar and the corresponding field will be denoted 
by 0. 

The Lagrangian describing the transition of a heavy quark Q into a lighter quark 
q and a "photon" has the form 

= hQ(j)q + h.c. , (2.9) 

where h is the coupling constant and Q = . The masses of the quarks Q and q are 
both large (and I remind that they are both spinless). Moreover, to further simplify 
the problem we will analyze a special limit (the so called small velocity or SV limit 
suggested in Ref. [|1^) in which 

Aqcd <^ rriQ - rUg TTiQ . (2.10) 

The field (j) carries color charge zero; the reaction Q ^ q + (j) could be considered a 
toy model for the radiative decays of the type B Xs'y where Xs is an arbitrary 
inclusive hadronic state containing the s quark produced in the b quark decay. 
It is very easy to calculate the total width for the free quark decay Q ^ q + (j), 

Tfree quark(Q 90) = ' T = Tq (2-11) 

SirrrLn 



where 



2 2 

E^ = —r> -■ (2-12) 



This free quark expression is valid for the total inclusive probability in the asymp- 
totic limit when mq —>■ oo. We are interested, however, in the preasymptotic cor- 
rections proportional to powers of l/mq. 

First of all we must formulate what object we must deal with in order to be able 
to calculate these corrections systematically. Upon reflection one concludes that it 
can not be the decay amplitude Q qcj) itself. Instead we must consider the Q Q 
forward "scattering" amplitude depicted on Fig. 1. By scattering I mean that Q 
scatters off the quantum and off the background gluon field which is not shown 
on Fig. 1 explicitly but is implied. It is implied that all quark lines, Q and g, are 
submerged into this soft-gluon background field. Through the optical theorem the 
imaginary part of the amplitude of Fig. 1 is related to the inclusive probability of 
the Q q(f) transition. More specifically, if we introduce the transition operator 

f = ij d^xe-''^''T{Q{x)q{x) , q{0)Q{0)}. (2.13) 
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then the energy spectrum of the particle in the inclusive decay is obtained from 
T in the following way: 



dV h^E , , " , . . . 

— = Im (Hn T Hn) . 2.14 



Here as usual Hq denotes a hadron built from the heavy quark Q and the light cloud 
(including the light antiquark), q in the exponent is the four-momentum carried away 
by (f) and E is the energy of the (p quantum. 

Equation ( p.l4| ) immediately translates the 1/mQ expansion for the transition 
operator in the l/rriQ expansion for the inclusive decay rate. The fact that the 
transition operator must be the primary object of the analysis in all problems of 
this type was realized in Refs. |P, 

Now we use what we have already learnt about the background field technique 
to write the transition operator in the form 

f = -Jd'xe-^''^{x\Q:^^Q\0) = 

- I dHAQji^ ^Q|0) , (2.15) 

where 

(Po)m = ruQVf, . 

The Green function of the quark q differs from the Green function given in Eq. ( p.6| ) 
in an obvious way since we assume for the time being that our quarks Q and q have 
spin zero, and, correspondingly, instead of Eq. ( [^.6D referring to the spinor quarks 
we have 

Six,0) = ix\ J O). (2.16) 

— rrig 

In the second hne of Eq. ( |2.15|) we proceeded to the rescaled fields Q which singles 
out the mechanical part of the momentum operator. 

One more thing which will be needed is the equation of motion for the scalar 
field Q substituting the Dirac equation. Starting from (V"^ — mQ)Q = we obviously 
get ^ 

^oQ = -^Q. (2.17) 

2mQ 

Finally we are ready to begin constructing the 1/mQ expansion. Since vr is 
of order Aqcd while Pq — q scales as itlq, in the leading approximation vr in the 
denominator of Eq. ( p. 15 ) can be neglected altogether. Then, obviously. 



= (2.18) 

mi — k"^ 
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where 

k = P^-q. 

We see that the leading operator appearing in the expansion is QQ and it has 
dimension 2 (let us recall that the scalar Q field has dimension 1 in contrast to the 
real quark fields of dimension 3/2, which leads in particular to different normalization 
factors in the matrix elements). Taking the imaginary part we conclude that in the 
leading approximation 

ilm {HQ\f^'^\HQ) = ^6iE - Eo) . (2.19) 
TT 2mQ 

Here and below I will use a very convenient short-hand notation, 

{QQ) ^ {Hq\QQ\Hq) . 

The delta function in the imaginary part is characteristic of a two-body decay. 
As a matter of fact, combining Eq. (|2.19| ) with the general expression (|2.14| ) and 



approximating (QQ) by unity - which can and must be done in the leading order in 
l/mg - we get the delta-function spectrum of the the free quark decay. Integrating 
over the energy we then arrive at the free quark decay width ( |2.11| ). 

Although this little achievement is quite gratifying and shows that we are on the 
right track the real l/mq expansion begins when the preasymptotic terms switch 
on. To this end the terms with vr in the denominator of Eq. ( |2.15|) must be kept. 



and then the expansion in [kn + vr^) must be carried out. The general term of this 
expansion is 

In Lecture 4 where the theory of the end point spectrum will be presented we will 
need the whole sum. At the moment our purpose is more limited - we are aimed at 
getting the first correction in the total decay width. This task does not require the 
infinite sum; only two terms, with n = 1 and n = 2, are relevant. Both terms are 
especially simple. 

Indeed, if n = 1 the combination 2mQ7TQ + vr^ in the numerator acting on Q is 
nothing else than the equation of motion, and can be dropped, see Eq. ( p.l7|) . We 
can further discard the gvr part - since the Hq spin is assumed to be zero there is 
no preferred orientation and, hence, {QttQ) = 0. In this way we arrive at 

{Hc,\f^'^\Hq) = f\ (2.21) 
rriQ [m^ — k'^y 



plus terms of higher order in l/mq. Here the same equation of motion ( p.l7|) was 
applied to eliminate ttq in favor of vr^ which, in the given order in l/mg, coincides 
with —7?^. The notation is even more concise than previously, namely (vr^) stands 
for {Hq\Q'k'^Q\Hq) . You will often see similar short-hand below. 
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I pause here to make a side remark. The physical meaning of the matrix element 
(tt^) is quite transparent - it merely represents the average value of the square of 
the momentum of the heavy quark Q inside the heavy hadron Hq. This quantity is 
of order Aqqj-, . This is one of the most important parameters of the heavy quark 
theory, along with A. Note the gap in dimensions of the operators appearing in the 
expansion. The dimension-2 operator QQ is followed by dimension-4 operator Qn^Q. 
No relevant operator of dimension 3 exists. Due to this reason the contribution of 
T^^^ in the total width is "unnaturally" suppressed by two powers of the inverse 
heavy quark mass, not one power as one would expect apriori. The observation of 
the dimension gap was first made in Ref. in the context of HQET; it is crucial 
in phenomenological applications. 

Let us return now to the construction of the l/mg expansion, and consider the 
term with n = 2 in the sum (|2.2(]| ). One of two factors (2mQ7ro + vr^) can be applied 



to the right, the other one to the left. The difference between vr applied to the 
right and to the left is a total derivative which vanishes anyway in the forward 
matrix element {Hq\...\Hq) . This simple observation implies that the combination 
(2mQ7ro + vr^) in the numerator again vanishes by virtue of the equation of motion 
and we are left with 

where I have singled out and retained only the spin-0 part of the operator QniiijQ 
{l/3)6ijQ7i^Q for the reasons explained above. 

We are almost done. The imaginary parts of T^^^ and T*^^) are expressible in 
terms of the first and second derivatives of the delta function, and after some simple 
algebra it is not difficult to get 

where operators of higher dimension are ignored; I have taken into account that 
qo = E and (f = E'^ and played a little with the delta functions. 

The expansion of ImT into local operators generates more and more singular 
terms at the point where the spectrum would be concentrated in the free quark 
approximation. You should not be surprised by this circumstance which will be 
elucidated in every detail in due time. What is important is that the physical 
spectrum is a smooth function of E. One could derive a smooth spectrum by 
summing up the infinite set of operators in Eq. ( |2.2CI| ) - this will be the subject of 
Lecture 4. There is no need to carry out this summation now, however, since we 
are interested only in the integral characteristics of the type of the total probability. 
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As far as such integral characteristics are concerned, the expansion in Eq. ( |2.23| ) is 
perfectly legitimate. 

At first, we calculate the total width by substituting eq. ( |2.23| ) into eq. ( P^.14D 
and integrating over E, 

where the integration runs from to the physical boundary Eq^^^ , expressed in terms 
of the hadron masses 

Ml -Ml 

° 2Mhq ^ ' 

The power correction proportional to (Qif^Q) / m?Q which might have appeared can- 
cels in the total width! Is this cancellation unexpected? No, we could have antic- 
ipated it on general grounds. Indeed, the total width F is a Lorentz scalar, and, 
quite naturally, the l/mq expansion for this quantity must run over the Lorentz 
scalar operators; QQ is Lorentz scalar while Qtt^Q is not. The fact that there are 
no explicit l/mg corrections in Eq. ( p.24| ) does not mean that they are absent in F 
at all. They could appear through {QQ)mQ/MHQ- Hence, our next task is to find 
the expansion for (QQ) in the toy model at hand. To solve the problem we will use 
the very same idea as in Sect. 1.5; the only difference is the form of the heavy quark 
current. For the scalar quarks the current whose diagonal matrix element counts 
the number of quarks is QiDfj.Q- Hence, in the rest frame of Hq we have: 

^ {HQlQiDoQlHg) = 1 . (2.26) 



Passing to the rescaled fields we arrive at the relation 

1 = j^{HQ\mQQQ + Qt^oQWq) = 
Mhq 

^{Hq\QQ\Hq) + — -i {HQ\Qn'Q\HQ), (2.27) 

where the second line is due to the equation of motion. Equation ( p.27| ) leads us to 
the result sought for, 

{Hq\QQ\Hq) =(i-St + a ; (2-28) 



I have introduced here the standard notation for the expectation value of vf^, 

„o.., = 5j^(^f«i2™Qe*=ei^'o> • (2-29) 
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As was already mentioned, /x^ is a crucial parameter of the heavy quark theory. 
Its definition in QCD will be slightly different from that in our toy model, but the 
physical meaning will be the same. 

Plugging in Eq. ( 2.28| ) in ( p.24 ) we finally arrive at the desired expression. 



The inclusive width coincides with the parton-model result up to terms of order 
l/mg. There is no correction of order l/mg! Moreover, the term l/rriQ is calcu- 
lable and its physical meaning is quite transparent: it refiects the time dilation for 
the moving quark inside the heavy hadron at rest (the Doppler effect). The coeffi- 
cient (—1/2) in front of (7f^)/mQ could, therefore, have been guessed from the very 
beginning, without explicit calculations, were we a little bit smarter. 

This situation is quite general, it takes place not only in the toy model at hand 
but in real QCD as well. The absence of the correction of order l/mg in the total 
inclusive widths (say, the semileptonic width of the B mesons, or T{B —>■ Xg'j), and 
so on) is called the CGG/BUV theorem fS^, |2§. 



I hasten to add, though, that the absence of the l/mq correction is not merely 
a consequence of the dimension gap in the set of the relevant operators, as it is 
sometimes stated in the literature. Indeed, let me give a counterexample. Let us 
calculate the average energy of the particle, or, more exactly, the first moment of 
the spectrum, 

= ,E(Et'--E)--. (2.31) 

In the parton model the spectrum is a pure delta function and, consequently, Ji 
vanishes. The heavy quark expansion does generate a non- vanishing result, a l/mg 
effect. To see that this is indeed the case we integrate the theoretical spectrum 
( ^.23| ) which yields us 

,,2 pphys 

/i = A - (2.32) 

where A is defined as follows 



- Eo, (2.33) 



and the parameters Eq and Eq'^^'^ are given in Eqs. (|]TD and (121251) . Now, invoking 



what we have already learnt in Sect. 1.5 about the heavy hadron masses we find 
that 

A = ^vlA + O(A^cd) (2.34) 

where 
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In the SV limit Vq is small and coincides with the velocity of the final hadron 
produced in the transition Q — > g0. 
With all these definitions 

/i = ^i;o'A + 0(A|cd) (2.35) 

i.e. the preasymptotic correction in the first moment is of the order of Aqcd? not 
^QCD- The reason for the occurrence of a "wrong" power of the QCD parameter is 
that the leading correction term in the 1 / mq expansion in this particular quantity 
is unrelated to any local operator. As we will see later on such a situation is not 
rare in the heavy quark theory. The sum rule ( p.35| ) is just a version of Voloshin's 
optical sum rule |3l|, while that of Eq. ( p.30|) can be interpreted in terms of the 



Bjorken sum rule [^. We will dwell on the both sum rules in the real QCD in Sect. 
3.6. 

I apologize for this little waterfall of new letters and definitions and hope that 
a simple picture behind our results is not overshadowed. Notice that in the SV 
limit Eq^^'^ — E reduces to the excitation energy of the final hadron produced in 
the decay. The factor Eq^^^ — E in the integrand eliminates the "elastic" peak, so 
that the integral is saturated only by the inelastic contributions. Say, in the 6 — > c0 
transition the contribution of -B — > D(j) is eliminated, only the excited D mesons 
survive in the first moment. Since the excitation energies are of order of Aqcd? ot A, 
the prediction ( f^.35| ) means that the probabilities of the inelastic transitions B 
excited D's are all proportional to v"^. This is in full agreement with the theorem 
T5| discussed in Sect. 1.3 - that in the point of zero recoil the only transition that 



can occur is the elastic B ^ D transition, with the unit probability. Away from the 
point of the zero recoil (but in the SV limit) the inelastic transitions are generated. 
However, Eq. ( |2.30|) shows, that up to small corrections 0(AQci-,/mQ) which can be 
neglected if we are interested only in the linear in Aqcd effects, the total probability 
remains unity. In other words, the total probability is just reshuffled: a small 
part is taken away from the elastic transition and is given to the inelastic transitions. 
The QCD analog of this assertion is the essence of the Bjorken sum rule [P^j. 

It is quite evident that the series of such sum rules can readily be continued 
further. For the next moment, for instance, we get 

Analyzing this sum rule in the SV limit one obtains, in principle, additional infor- 



mation, not included in the results of Refs. ||T5|, |3T|, |3^ . It is worth emphasizing that 
in Eqs. ( |2.24D , ( |2.32D and ( ^^.361 ) we have collected all terms through order Aq^j-,, 
whereas those of order Aqqj-, are systematically omitted. Predictions for higher 
moments would require calculating terms 0{Aqqy)) and higher. 

Concluding this part let me suggest to you an exercise which will show whether 
the technology introduced above is well understood by you. Try to repeat in real 
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QCD, with the quark spins switched on, everything we have done in the toy modeL 
Of particular interest to us will be the transition operator 



T 



d^xe-''i^T{Q{x)T^q{x) , g(0)r,Q(0)}. 



(2.37) 



where is either 7^ or 7/^75. This transition operator is relevant for the semileptonic 
6 to c decays. To facilitate the task consider special kinematics: (i) zero recoil (the 
vanishing spatial momentum of the lepton pair, g* = 0; (ii) small velocity limit 
(f <^ mq. To further facilitate the task limit yourself to the spatial components of 
. If you still have problems go over this lecture again and consult the original 
55| , The full answer for the transition operator ( |2.37D is given, for 

instance, in Appendix of Ref. . 



works 



2.2 The Fock-Schwinger gauge 

In some situations (especially when one deals with massless quarks) a variant of 
the background field technique based on the so called Fock-Schwinger gauge for the 
external filed turns out to very efficient (for a review and extensive list of references 
see |^8[). The gauge condition on the background gluon field has the form 

XpAp{x) = 0. (2.38) 

What is remarkable in this condition is that in this gauge the gauge four-potential 
can be represented as an expansion which runs only over the gauge covariant quan- 
tities, the gluon field strength tensor and its covariant derivatives, 

Ap{x) = ^XpGp^iO) + i^x«XpD„Gp^,(0) + - • (2.39) 

This expression implies, in particular, that ^4(0) = 0. It is worth noting that the 
gauge condition ( |2.38| ) singles out the origin and, hence, breaks the translational 
invariance. The latter is restored only in the final answer for the gauge invariant 
amplitudes. 



It is rather easy to show (see p8|) that the massless quark Green function in the 
coordinate space is 

'S'(x, 0) = - Ga07075 + Ga<p = -ea,i>fj.uGf,iy. (2.40) 

One can also construct a similar expansion for the Green function in the momentum 
space S{q). 

If the quark is not massless, rriq 7^ 0, the expansion of the Green function in the 
background field becomes much more cumbersome. Although we will hardly need 
it in full I quote it here for the sake of completeness. 



1 y 1 2 2j^/ / 2\\ 1 ^"/=< j -x^mKi{mV^^ 

S{x,0) = TTT^j-TT"^ ^ K2{mV-x^)} - ^^Ga07075 



271^ x^' 2 ' 8n^x^ 
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^ r\ -I I / / / A / I — I 1 1 1 

pxcj' i^o(W-x2) + .... (2.41) 



167r2 



Here K is the McDonald function. This result was obtained in Ref. Further 
education on the Fock-Schwinger gauge technique can be obtained from Ref. p8| . 
The best way to master those aspects which are most common in the heavy quark 
theory is merely to play with this tool kit. Let us see how it works in the calculation 
of the l/rriQ correction in the total probability of the semileptonic decay of the heavy 
quark in the real QCD. Unlike Sect. 2.1 we will address directly the total width 
bypassing the stage of the spectrum. The first calculation of the power correction 
in T{B XJu) along these lines was carried out in Ref. ||23 (see also p7[]). 



2.3 The l/ruQ corrections to the semileptonic inclusive width 
in QCD 

In this section I will describe probably the most elegant application of the ideas 
developed above - calculation of the leading correction in the total semileptonic 
widths. In the toy model considered in the previous section it was established 
that the 1/mQ correction was absent, and the first non-trivial correction l/rriQ was 
associated with the matrix element of (QQ)- As a matter of fact, through the heavy 
quark expansion, we managed to express it in terms of the matrix element (Qif'^Q). 

When we take the quark spins into account a new dimension-5 Lorentz scalar 
operator appears, {Q{i/2)aGQ) . On general grounds one may expect that the main 
lesson abstracted from the toy model - the absence of the l/ruQ term - persists 
but the ^/ttlq correction will receive a contribution from the operator Q{i/2)aGQ. 
The conclusion will be confirmed by the analysis presented below. You will see how 
efficiently the Fock-Schwinger technique is in this case. 

Thus, let us proceed to calculation of semileptonic widths. The final quark mass 
ruq is arbitrary - we do not assume the SV limit now, nor is any other constraint 
imposed on rrig. The weak Lagrangian responsible for the semileptonic decays has 
the following generic form 

C = ^VQ,{qT,Q)ilT,u), r^ = 7^(l+75), (2.42) 

where / is a charged lepton, electron for definiteness. The mass of the charged lepton 
will be neglected. Moreover, Gp and Vqq are constants irrelevant for our purposes. 
As usual, at the first stage we construct the transition operator T{Q X —> Q), 

f = ij d^xT{C{x)C{<d)} = CiOi (2.43) 

describing a diagonal amplitude with the heavy quark Q in the initial and final 
state (with identical momenta). The lowest-dimension operator in the expansion of 
T{Q X Q) is QQ, and the complete perturbative prediction - the spectator 
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model - corresponds to the perturbative calculation of the coefficient of this operator. 
For the time being we are not interested in perturbative calculations. Our task is the 
analysis of the influence of the soft modes in the gluon field manifesting themselves 
as a series of higher-dimension operators in T. 

At the second stage we average T over the hadronic state of interest, say, B 
mesons. At this stage the non-perturbative large distance dynamics enters through 
matrix elements of the operators of dimension 5 and higher. 

Finally, the imaginary part of {Hq\T\Hq) presents the Hq semileptonic width 
sought for, 

T = -^lm{HQ\T\HQ) . (2.44) 

Mhq 

The diagram determining the transition operator is depicted on Fig. 2. The 
lepton propagators are, of course, free - they do not feel the background gluon field. 
Thick lines refer to the initial quark Q. Although the gluon field is not shown one 
should understand that the lines corresponding to Q and q are submerged into a 
soft gluon background. 

In the Fock-Schwinger gauge the line corresponding to the final quark q (Fig. 2) 
remains free, and the only source of the dimension-5 operators is the external line 
corresponding to Q (or Q). Let us elaborate this point in more detail. 

If we do not target corrections higher than l/mg it is sufficient to use the ex- 
pression for the quark Green function given in Eq. ( p. 401 ) or ( |2.41| ). The particular 



form is absolutely inessential; the only important point is the chiral structure of the 
vertices in the weak Lagrangian and the fact that the leptons are massless. 

The currents in the weak Lagrangian ( p.42|) are left-handed. Therefore, the 



Green function of the quark q is sandwiched between and F,^. This means that 
1 + 75 projectors annihilate the part of the Green function with the even number 
of the 7 matrices. Then the only potential contribution is associated with the ffist 
line in eq. ( |2.41| ). 

The non-perturbative term is the one containing Gaif,- This term vanishes, how- 
ever, after convoluting it with the lepton part. Indeed, the lepton loop (with the 
massless leptons) has the form 

2 1 

L^,u = j^{2x^x^ - X g^^). (2.45) 

Here I take the product of two massless fermion propagators in the coordinate space, 
with the appropriate 7 matrices inserted, and do the trace. Actually we need to 
know only the last bracket. Now, convoluting it with the G term from the quark 
Green function we get 

{T^XaGa<j>'l^Ty){2x^Xu - X^Q^y) = 0, 

q.e.d. 0,0. 
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Thus, if one uses the Fock-Schwinger gauge the only source for the l/mq cor- 
rection in the total semileptonic widths (at the level up to I/ttlq) is through the 
equations of motion for Q. Here is how it works. 

The expression for the amplitude corresponding to the diagram of Fig. 2 can be 
generically written as follows: 

A = J d^xQ{0)F{x)Q{x), (2.46) 

where a function F{x) incorporates the lepton loop and the q quark Green function. 
It may include Lorentz and color matrices, etc. Now, let us single out the large, 
mechanical part of the motion of the heavy quark, 

Q(x) = e-^^o^g(x). 

Then 

A = Jd'x |0(0)F(x)e-^^o^[Q(0) + x^d^Q{0) + ^x^x^d^dM^) + ■■•]} = 
|(?(0)[F(Po)Q(0) +^^F(Po)5mQ(0) 

0(0)F(Po + id)QiO) = Q(0)F(z9)Q(0), (2.47) 

where F is the Fourier transform of F{x). 

Our next goal is to convert id in the covariant derivative and then use the 
equation of motion, i PQ = rnqQ. More exactly, we start from the expressions of 
the form 

Q(0) Ap'fQi.^). P, = ^d^. (2.48) 

rewrite in terms of = iD^ plus terms with the gluon field strength tensor (in 
the Fock-Schwinger gauge) and then substitute 1/ acting on Q by mq. Expressions 
( p.48| ) appear in ImT. If the final q quark is massless, niq = 0, the only relevant 
power is /c = 2. Switching on the quark mass, ruq ^ 0, brings in other values of k 



as well. (Warning: in the procedure sketched above all operators p in Eq. ( p. 48 
should be considered as acting either only to the right or only to the left. I will 
assume they act to the right. We can not make some of them act to the right and 
others to the left and neglect full derivatives. Question: do you understand why?) 

Since we focus now on QaGQ it is sufficient to keep only the terms linear in the 
gluon field strength tensor; the terms with derivatives of G^i, are to be neglected as 
well. In this approximation in the Fock-Schwinger gauge = {l/2)xpGpa- 

Furthermore, 

/ = p2 - 2Ap (2.49) 
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where we neglected the terms quadratic in A and used the fact that [p^, A^] = 0. 
Equation (|2.49|) should be substituted in Eq. ( |2.48| ), and then we start transposing 
Ap trying to put it in the left-most position, next to Q(0). If Ap appears in this 
position the result is zero since ^4(0) = 0. Notice that 

[p^, Ap] OC Ga/SPaPlB = 0, 

SO, one can freely transpose Ap through p'^. In this way we arrive at 

Q(0) Mp')'Q(O) = Q(0) [nv')" - 2k f{AV){V^f-^] g(0) . (2.50) 
Moreover, with our accuracy the last term reduces to 

= -l{f<yG -aGf)f^^-^ (2.51) 
8 

and, hence, using the equations of motion we conclude that 

Q(0) f{AV){r^f-^Q{'d) = 0. (2.52) 
As a result, the aG terms emerge only due to the fact that 

v'=r-'^cTG, 

and the final expression is as follows: 

Q(0) p[p')'Q{0) - m5+iQ(0)Q(0) - jQiO)aGQ{OWQ-\ (2.53) 

After these explanatory remarks the procedure of calculating the leading l/rriQ 
correction in the total semileptonic width should be perfectly clear. Let us summa- 
rize it in the form of a prescription. 

(i) Calculate the semileptonic width in the parton model. The result has the 
form 

^-4^0(5) (2.54) 

where Fo(mg/mQ) is the phase space factor, a function of the ratio m^/mg well- 
known in the literature (see Eq. ( |2.57|) below). 

(ii) Then construct the expression for T including the 0{'m~^) corrections in the 



following way |24 
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+ 4 iP^ - m{p) \ (2.55) 



where 



= ,^{HQ\Qit/2)aGQ\HQ) (2.56) 
2M/^Q 



and 

P 



A few comments are in order concerning this beautiful expression for the total 
semileptonic width. The expansion contains two Lorentz scalar operators, QQ and 
Q{i/2)aGQ, of dimension 3 and 5, respectively. The fact that only the Lorentz 
scalars contribute is obvious since F is a Lorentz invariant quantity. We observe 
here the very same gap in dimensions mentioned previously in the context of the 
toy model - there is no operator of dimension 4 [^]. The only element still needed to 
complete the derivation is the matrix element (QQ) . Fortunately, the corresponding 
heavy quark expansion has been already built, see Eq. ( |1.3CI| ). 

Borrowing the explicit expression for -Fo(p) from textbooks (it is singled out 
below in the braces) and assembling all other pieces together we finally get 

192 7r3 ^ 



where /i^ is defined in QCD as follows 

pI = 77^{Hq\Q^'Q\Hq). (2.58) 

This result is due to Bigi et al. P3]. The absence of the l/rriQ correction is a 
manifestation of the CGG/BUV theorem. 

If the mass of the final charged lepton is non- negligible, the property of no 
soft gluon emission from the q quark line (in the Fock-Schwinger gauge) is lost. 
The expansion for F(i? XcTUt) in this case is much more cumbersome; it was 
constructed in Ref. . 




(2.57) 



Dimension-5 operators are responsible for the leading non-perturbative correc- 
tions in the total semileptonic widths. To assess the convergence of the expansion it 
may be instructive to have an idea of the higher order terms in the expansion. The 
I/ttlq terms were estimated in Ref. |^D[. This is a rather messy and time-consuming 



analysis, and it is hardly in order to comment on it in this lecture. Surprising though 
it is, from what we already know it is practically trivial to find the coefficient of the 
dimension-7 operator 

QG^.G^.g, (2.59) 
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with two gluon field strength tensors fully contacted over the Lorentz indices for 
massless final quarks (say, b ^ u transition). Of course, this is a purely academic 
exercise, for many reasons. In particular, because it is only one of a rather large 
number of dimension-7 operators. Since we do not know their matrix elements 
anyhow, it seems to be meaningless to carry out full classification and calculate all 
coefficients. The operator (|2.59| ) is chosen since one can at least use factorization for 
a rough estimate of the corresponding matrix element and since we get its coefficient 
essentially for free. 

The point is that the massless quark propagator (in the Fock-Scwinger gauge) 
does not contain GuyGfj^^ term at all (see ref. [0). This fact implies that the only 
source of the operator (|2.59|) is the same as in the case of QcrGQ. It is not difficult 



to get that 

Q(o)#g(o) = g(o) fv^Qio) + ^Q(o)G^,G^,g(o) = 

Q{0)f'Q{0) + Q{0)G,,G,MO), (2.60) 
where we omitted structures of the type [G^j^^G^a — {^/4:)guaG^pG^p]. 



2.4 Digression 

This section is intended for curious readers - those who are anxious to find out where 
and how else, beyond the theory of the Hq states, the background field technique 
can be used to obtain interesting predictions. Here I will discuss an estimate of the 
mass splittings between the levels of the highly excited quarkonium states. This 
part can be safely omitted in first reading since it is unrelated to the remainder of 
these lectures. 

The quarkonium states to be considered below consist of one quark Q, one anti- 
quark Q and the soft gluon cloud connecting them together. To begin with we will 
assume that mq is large, niQ ^ Aqcd (later on this assumption will be relaxed). 
The QQ mesons can have different quantum numbers. We will analyze the excited 
S and P wave states with the quantum numbers 0~ and 0^, respectively. The naive 
quark model language is used to name the states; this does not mean, of course, that 
we accept any of the dynamical assumptions of the naive quark model. It would be 
more accurate to say that the mesons of interest are produced from the vacuum by 
the currents 

Jp = Qi-i^Q and Js = QQ ■ (2.61) 

The central object of our analysis is the difference between the two-point func- 
tions in the pseudoscalar and scalar channels. In terms of the Green functions in 
the background field this difference takes the form (see Fig. 3) 

Up-Us = i [ dxe''^'' {{vac\T{Jp{x)Jp{0)}\vac) - {vac\T{Jsix)JsiO)}\vac)) = 
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1 1 1 



A comment is in order here concerning the trace operation in this expression. It 
imphes not only the trace over the Lorentz and color indices, as usual, but also the 
trace in the momentum space substituting / rf^p/ (27r)^ in the conventional Feynman 
integral. With the help of Eq. ( p.6|) the difference lip — fl^ can be identically 
rewritten as 

Hp - = -2m^zTr [^^^ (2-63) 

where 

D^ = [V±{q/2)f~ml + {i/2)aG. 

We continue to ignore hard gluons assuming that the only role of the gluon field 
is to provide a soft cementing background. This is certainly an idealization, but let 
us see how far one can go within the framework of this simplified picture. Neglecting 
hard gluons means, in particular, that we will be unable to analyze the low-lying 
levels of heavy quarkonium where an essential role is played by the short- distance 
Coulomb interaction. Each aG insertion in the denominator is of order of Aqqj-,, i.e. 
does not scale with the external momentum q when q is large. Let us expand Eq. 
( p.63| ) in aG and take the trace over the Lorentz indices. Then the first order term 
drops out; the first surviving term is bilinear in G, 

Hp - = -8m|iTr i ^ ^ + 



11 1 1 111^ 1^ 1 1 

+ (2.64) 

A closer look at this expression reveals some peculiar features. First of all, one 
can interpret each term as a certain correlation function in the theory where the 
quark Q is scalar, not spinor. Take, for instance, this first line. It is nothing else 
than the two-point function of the L = quarkonium in the scalar QCD (i.e. QCD 
with the scalar quarks; L is the total angular momentum of the meson). The current 
producing the scalar quarkonium from the vacuum in the scalar QCD is Q^Q (Fig. 
4). The second and the third line, together, represent the four-point function of the 
type depicted on Fig. 5. The current denoted by the dashed line on this figure is 
Q^GaisQ] the momentum flowing through this line vanishes. Two insertions of G 
imply that this four-point function is explicitly proportional to Aq^^q. 

Now, let us examine Eq. ( |2.64| ) in the complex q"^ plane. At some positive values 
of the two-point function of Fig. 4 has simple poles corresponding to positions of 



32 



the L = quarkonium levels in the scalar QCD. The four-point function of Fig. 5 
has double and single poles at the very same values of and single poles at some 
other values of corresponding to the production of L = 1 states in the scalar 
QCD. The latter are produced from the L — states by applying to them the 

current Q^GapQ- 

On the other hand, the original difference lip — in real QCD has only single 
poles at the positions of the S and P wave states. These positions are shifted 
compared to the levels in the scalar QCD. Expanding in the shift one generates 
double poles. The L — 1 pole - the partner to the P wave meson states in lip — II^ 
appears only in the four-point function of Fig. 5. From this figure it is quite clear 
that the residue of the L = 1 pole scales as Aqqj-,/(A(M^))^ where A(M^) is a 
characteristic splitting between the L = and L = 1 states. On the other hand, 
the residue of the P wave meson in lip — n^, on general grounds, scales as p^/M^ 
where p is a characteristic quark momentum, and I assume that Aqcd <S |p| <S M. 
Equating these two estimates we find that 

AM~A^cd/Ip1- 

One may observe, with satisfaction, that this is exactly the characteristic level split- 
ting (between radial or orbital excitations) for two heavy quarks interacting through 
a string ("linear potential"). What is remarkable is that in no place our estimate 
invokes any reference to the linear potential or other models. It was based only on 
some general features of QCD. For me this is a strong evidence that a string-like 
picture should take place in QCD, at least, approximately, for high excitations. 

What changes if the quarks are light or even massless, itlq oo? The only 
difference is that now the residues of the P wave mesons in lip — II^ are of the same 
order as those of the S wave mesons for highly excited states, which implies that 

A^qcd/(A(M2))2 ~ 1 , 

or 

A(M2) ~ A|cD ■ 

In other words, we got the linear Regge trajectories, at least for highly excited states. 
Moreover, this analysis makes clear a potentially important point - the empirical 
observation that even the lowest states in every channel lie on the linear Regge 
trajectories looks like a numerical coincidence and can not be exact. 
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3 Lecture 3. Classic Problems with Heavy Quarks 



The number of problems successfully solved within the heavy quark expansion is 
quite large. Even a brief review of the main applications is beyond the scope of these 
lectures. Some issues, however, are quite general and are important in a variety of 
applications. Everybody, not only the heavy quark practitioners, should know them. 
In this lecture we will discuss several such topics - the scaling of the heavy meson 
coupling constants, some properties of the Isgur-Wise function, and, finally, analysis 
of corrections violating the heavy quark symmetry at the point of zero recoil. We 
begin, however, from a systematic classification of all local operators which appear 
in the heavy quark expansion up to the level 0{mQ^). Some terms of order l/m^ 
in particular heavy quark expansions are actually not expressible in terms of the 
local operators and are, rather, related to non-local correlation functions. A full 



classification of such correction also exists |42, 25], but we will not go into details 



only marginally mentioning them here and there. The interested reader is referred 
to the original publications B3]. 



3.1 Catalogue of relevant operators 

The local operators in the heavy quark expansion are bilinear in the heavy quark 
field. They are certainly gauge invariant, and in many instances, when the expan- 
sion is built for scalar quantities, the operators must be Lorentz scalars. As in any 
operator product expansion in QCD they can be ordered according to their dimen- 
sion. We will limit ourselves here to dimension 6 and lower. This leaves us with 
quite a few possibilities listed below. We start with the Lorentz scalar operators. 
The only appropriate operators are 

_ 1 _ _ 

QQ , -Qaf^yG^yQ , and QVQqVq (3.1) 

where F stands here for a combination of 7 and color matrices. All other structures 
that might come to one's mind reduce to those listed above and full derivatives by 
virtue of the equations of motion. 

(Exercise: prove that this is the case, for instance, for the operators QD^Q and 

(i) The only operator of dimension 3 is QQ. This operator is related to the heavy 
quark current Q'joQ plus terms suppressed by powers of l/mg. The leading term of 
this expansion has been already discussed, see Eq. ( |1.3CI| ). Actually it is not difficult 
to continue this expansion one step further. The following relation is exact: 

QQ = QioQ + 2gfi^)'g = Q7og - 2Q-^--^q = 

\ 2 / Inin 2mn 



QI0Q + QtzttQ + ^ total derivative ; (3.2) 



2m: 



Q 
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in the first relation above the operators vr^ act on the Q field. Keeping in mind 
that we always consider only the forward matrix elements, with the zero momentum 
transfer, we can drop all terms with total derivatives. Applying now the equations 
of motion ( |1.8|) and (|1.9|) generates the l/wig expansion for the scalar density. 



(3.3) 



Here = GiQ is the chromoelectric field, and its covariant derivative is defined as [] 
DjEk = —i[7ij, Ek]; we have omitted the term Q{[nk, [ttq, VTj]] — [vTj, [vro,7rfc]])Q using 
the Jacobi identity. Moreover, 

DE = gh^j;; 

by virtue of the QCD equation of motion (here J° = J2qQ1^it"'Q is the color quark 
current). Therefore the first of the l/rriQ terms can be rewritten as a four-fermion 
operator. 

(ii) As has been already mentioned, no operators of dimension 4 exist. 

(iii) Dimension five. There is only one such operator, 

Og = Q'-cr^uG^,Q, (3.4) 
where a^u = (l/2)[7|j, 7!^]- Again, it can be expanded in the powers of l/mq, 



Og = -QaBQ - -^Q (-{DE) + 2a ■ E x n) Q + O ] . (3.5) 

where B is the chromomagnetic field, B = V x A = —(l/2)eijkGjk- 

(iv) Dimension 6 four-quark operators J2i Q^iQQ^iQ- Generally speaking, the 
matrix Fj can be any Lorentz matrix (1,75,7^,7^75 or a^^) or any of the above 
multiplied by t". Of course, in specific problems only a subset of these matrices may 
appear. The four-quark operators differ by the chiral properties of the light quark 
field q. Some of them carry non-zero chirality (they are non-singlet with respect to 
SU{Nf)L X SU{Nf)fi). Hence, they do not show up in the transitions associated 
with the weak currents of the V — A type. 

Further remarks will concern operators that are spatial scalars but not Lorentz 



scalars. They appear in the low-energy effective Lagrangian (|1 . 1 1| ) and in the ex- 
pansions of the type (|3.3| ) and ( p.5| ). The most important operator from the class 
is 

a = Q^^Q , (3.6) 

which we have already encountered more than once. 



^Note that in our notations D = —d/dx — lA, therefore {DE) = — divE in the Abehan case. 
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Dimension 4 operators are all reducible to those of dimension 5 and higher. For 
instance, 

At the level of dimension 6 only one additional operator emerges (apart from the 
four-fermion operators), namely, 

Qa-ExnQ. (3.7) 

At first sight it might seem that one could build extra operators of dimension 6, 
from the gluon fields, e.g. 

QiTiEiQ or Qaitijk{DjEk)Q. 

Actually they are reducible to operators of higher dimension via the equations of 
motion. Indeed, using the fact that 

Ei = -ifTToTTj] 

one can rewrite 

QiTiEiQ = -iQiril-KoTTilQ = -iQ{'Ki'KQ'Ki - -Kf-K^jQ = 

11 

-i(5([7ri7ro]7ri + ttott,^ - 7r,^7ro)(5 ^ -iQ([7ri7ro]7ri + (o^vf) V,^ - 'K'^{anf)Q 

% — 

~i:Q[^i[^^'^i[]Q + dimension seven — 

— —iQ{diyE)Q + dimension seven. (3.8) 
This is a four-fermion operator. By the same token, 

Q(^i^ijk(DjEk)Q — QaiDoBiQ — 

- iQailTToBilQ = -7r^g[(a7f)2, ^^|g ^3 9^, 

2mQ 

which is obviously of the next order in I/uiq (a dimension-seven operator). 
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3.2 Extracting/determining the matrix elements 

Construction of the operator product expansion is only the first step in any theo- 
retical analysis. The heavy quark expansion must be converted into predictions for 
the physical quantities. To this end it is necessary to take the matrix elements of 
the operators involved in the expansion. The latter carry all information about the 
large distance dynamics responsible for the hadronic structure, in all its peculiarity. 
These matrix elements in our QCD-based approach play the same role as the wave 
functions in the non-relativistic quark models. 

In this section we will summarize what is known about the matrix elements of 
the operators from the list presented above. 

(i) The most favorable situation takes place at the level of dimension three. 
Indeed, the only Lorentz scalar operator of dimension 3 is QQ which has a nice 
expansion ( |3.3|) . The operator Q'joQ is the time component of the conserved current, 
measuring the number of the quarks Q in Hq. Therefore, both for mesons and 
baryons 

^ {Hq\Q^oQ\Hq) = 1. (3.10) 



2Mh 

(As usual, we stick to the rest frame of Hq; in the case of baryons averaging over 
the baryon spin is implied). 

(ii) The status of two operators of dimension 5 is different. Let us consider first 
Og whose matrix elements are expressible in terms of experimentally measurable 
quantities. 

To leading order in l/mq the parameter /i^ defined in Eq. ( ^.5(j| ) reduces to f\ 

/"G = 77^{Hq\Og\Hq) = --^{HQ\QaBQ\HQ) . (3.11) 

For pseudoscalar mesons this quantity can be related to the measured hyperfine 
mass splittings. Indeed, QaBQ is the leading spin-dependent operator in the heavy 
quark Hamiltonian (|1.14|) . Hence 

/x^(Sq) = ^(M|.-M|) (3.12) 

where B* and B are generic notations for the vector and pseudoscalar mesons, 
respectively, and the limit mq — oo is implied. Assuming that the h quark already 
belongs to this asymptotic limit one estimates ixq from the measured B meson 
masses, 

/i^ ^ 0.35 GeV^ . 

Furthermore, in the baryon family four baryons are expected to decay weakly 
and are, thus, long-living states: Ag, Eg, Sq and VLq. In the first three of them 

® Some authors prefer a different nomenclature p3| . The expectation values of the chromomag- 
netic and kinetic energy operators are sometimes called A2 and Ai. 
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the total angular momentum of the light cloud is zero; hence the chromomagnetic 
field has no vector to be aligned with, and 



/i^(AQ) = ^^(Sq) = filing) = 0. (3.13) 
In the case of Qq the total angular momentum of the light cloud is 1. Hence, 



l,U^Q) = ^[M'^,^,-M'^,^,j (3.14) 

where the superscripts 3/2 and 1/2 mark the spin of the baryon. Although the 
mass splitting on the right-hand side of Eq. (|3.14|) is in principle measurable, it 
is not known at present, and if one wants to get an estimate one has to resort to 
quark models or lattice calculations. Both approaches are not mature enough at 
the moment to give reliable predictions for this quantity and I suggest we wait until 
experimental measurements appear. 

Let us proceed now to the discussion of the matrix element of the operator 0.,^. 
The physical meaning of this matrix element is the average kinetic energy (more 
exactly, the spatial momentum squared) of the heavy quark Q inside Hq. This 
operator is spin- independent, and it is much harder to extract /i^ (the parameter 
fi^ is defined in Eq. ( |2.58| )) from phenomenology, although such an extraction is 
possible, in principle (see Ref. ||2^ and Lectures 4 and 5 for details). Since the 
phenomenological analysis has not been carried out yet one has to rely on theoretical 
estimates. Several calculations of fi^ within the QCD sum rules yield |4^, ^ 



l^liB) = ^(B\0^\B) = 0.5 ± 0.1 GeV^ . 

A remarkable model-independent lower bound on ^j^{B) exists in the literature, 

nl{B) > nl{B) ^ 0.35 GeV^ (3.15) 

The quantum- mechanical derivation of this inequality due to Voloshin (see Ref. 0]) 
is straightforward. Indeed, start from the square of the Hermitean operator (avf)^ 
and average it over the B meson state. It is obvious then that (((Ttt)^) > 0. Using the 
fact that (crvr)^ = vr^ + aB we immediately arrive at Eq. ( |3.15| ). A field-theoretic 



derivation of the same result can be found in Ref. p5[. It is remarkable that 
the inequality (|3.15| ) almost saturates the QCD sum rule estimate quoted above. 



Another lower bound on /i^(i?), obtained from a totally different line of reasoning, 
is discussed in Sect. 3.6. It turns out to be close to Eq. ( p.l5| ) numerically. 

It is plausible that for mesons and baryons is different - there is no reason 
why they should coincide. The task of estimating fi^. for baryons remains open. 

These parameters, fi^. and fiQ, along with A, are most important in applications. 
In most applications one deals with the expectation values over the B meson state. 
Therefore, let us agree that /i^, fi^, and A, with no subscripts or arguments, are 
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defined with respect to the B mesons. This is a standard convention. In a few cases 
when these quantities are defined with respect to other heavy flavor hadrons we will 
mark them by the corresponding subscripts or indicate with parentheses. 

(iii) Operators of dimension 6 are studied to a much lesser extent than those 
of dimension 5. Perhaps, the least favorable is the situation with the operator Oe 
given in Eq. ( |3.7| ). Let us parametrize its matrix element as follows: 

^ {HQlQa X EttQIHq) = 11%. (3.16) 



2Mh 

This operator in the heavy quark Hamiltonian is responsible for the spin-orbit in- 
teractions and consequently generates the spin-orbit splittings between the masses 
of the ground states and the orbital excitations. Hence, in the non-relativistic limit 
(non-relativistic with respect to the spectator light quark) /i^ vanishes for the S 
wave states. Of course, the non-relativistic approximation with respect to the light 
quark is very bad. The estimate of existing in the literature |0 is so rough that 
it, probably, does not deserve to be discussed here. 

As for the four-quark operators the only method of estimating their matrix ele- 
ments which does not rely heavily on the most primitive (and hence totally unreli- 
able) quark models is the old idea of factorization applicable only in mesons but - 
alas - not in baryons. 

First of all let us observe that each of the four-quark operators exists in two 
variants differing by the color flow. One can always rearrange the operators, using 
the Fierz identities, in the form 

Oig = QrqqrQ (3.17) 

and 

dig = Qt''Tqqt''TQ . (3.18) 

Take for definiteness F = 7^75. (Other 7 matrices can also appear, of course.) In the 
first operator color is transferred from the initial heavy to the final light quark and 
from the initial light to the final heavy quark. The second operator is essentially 
color-exchange. Now, if we are interested in the matrix elements over the meson 
states we can simply factorize the currents appearing in Eqs. (|3.17|) and ( p.l8| ) (i.e. 
saturate by the vacuum intermediate state), 

^{BQ\0,g\BQ) = ^MBfl ^{BQ\d,g\BQ) = (3.19) 

where Jb is the pseudoscalar decay constant, 

{0\Ql,lrMBQ) = ifBP^ (3.20) 

As we will discuss shortly, in the limit rriQ 00 the combination MbJb scales as 
ttiq (modulo logarithmic corrections) so that the right-hand side of Eq. (|3.19| ) is the 
cube of a typical hadronic mass, as it should be. 
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Factorization in Eq. ( p.l9| ) is justified by l/N^ arguments. Indeed, corrections 
to Eq. ( p.l9| ) are of tlie order of A*""^. 

Thus, from the whole set of the four-quark operators we can say something about 
the meson expectation values of those operators which are reducible to 

{QVq){qVQ) (3.21) 

where F stands for a Lorentz matrix but not for the color one, and the color indices 
are contracted within each of two brackets separately. Up to terms 0{N~'^) two 
brackets can be factorized. 

To get an idea about the numerical value of (1/2)/|Mb we should substi- 
tute a numerical value for Jb which is not measured so far. Theoretical ideas 
about this fundamental constant will be discussed later. Now let me say only that 
(1/2)/|Mb ~ 0.1 GeV^ with a significant uncertainty. Those matrix elements 
which are due to nonfactorizable contributions (see Eq. ( ^.181 )) are essentially unde- 
termined, although they are expected to be suppressed compared to the factorizable 
matrix elements (l/2)/gMs. 

As for the baryon matrix elements of the four-quark operators next-to-nothing 
is known about them at the moment. Some very crude estimates within the naive 



quark model are available but they are very unreliable. 

In conclusion of this section a remark is in order concerning numerical estimates 
of the key parameter of the heavy quark theory, A. I postponed discussing the issue 



because its value continues to be controversial. QCD sum rules indicate ^ 
that A ~ 0.5 GeV. This number is in full agreement with the lower bound stemming 
from Voloshin's sum rule, see Sect. 3.6. However, some lattice calculations yield a 
factor of 2 lower estimate. I am inclined to think that there is something wrong in 
the lattice results. Perhaps, the lattice definition of A does not fully correspond to 
that of the continuum theory. It is inconceivable that such a low value of A as 0.2 or 
even 0.3 GeV could be reconciled with the lower bound implied by Voloshin's sum 
rule. 

In the discussion above we have totally disregarded logarithmic dependence of 
the operators and their matrix elements due to anomalous dimensions - i.e. the 
issue of the normalization point (including the normalization point of A). This is 
in line with that so far we pretend that hard gluons do not exist. A brief excursion 
into this topic will be undertaken later; here it is only worth mentioning that all 
numerical estimates presented above refer to a low normalization point, of order of 
several units times Aqcd- 

We are ready now to review classic problems of the heavy quark theory. We will 
gradually move from simpler to more sophisticated problems. 

3.3 Mass formula revisited 

In Sect. 1.5 we have found the first subleading term in the mass formula for the 
heavy flavor hadrons. The parameter A was related to the expectation value of 
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the gluon anomaly, see Eq. ( |1.32| ). It is very easy to continue the expansion one 



step further and find the next subleading term, of order l/mq. One could have 
extended the derivation along the lines suggested in Sect. 1.5. This was done in Ref. 



25|. This is not the fastest route, however. Instead, let us observe that the l/mq 
term in the Hamiltonian ( |1.14| ) can be considered as the first order perturbation; 
the corresponding correction to the mass is merely the expectation value of this 
perturbation, 

MH^=mQ + A + ^{2Mh^)-'{Hq\^^ + aB\HQ) + ... = 

mQ + A + ^^^V^ + - (3-22) 

The terms of order I/tOq are neglected. If we keep only the terms up to l/mq 
it does not matter whether the state Hq we average over is an asymptotic state 
(corresponding to mq = oo) or the real physical heavy fiavor state. I remind that, 
unlike HQET, we work with the physical states. The difference becomes noticeable 
only at the level l/rriq. In this order the mass formula does not reduce any more to 
the expectation values of local operators. A part of the l/mq correction is due to 
non-local correlation functions, see Ref. for further details. Eq. ( |3.22[ ) was first 



presented in Ref. ES 



3.4 The scaling law of the pseudoscalar and vector coupling 
constants 

The pseudoscalar and vector meson constants fp and fv are defined as 

m-f.j^qlBq) = ifpp, mi,q\B*q) = ifyMye^ . (3.23) 

An alternative definition of the pseudoscalar constant can be given in terms of the 
pseudoscalar current, 

{Q\Qi^^q\Bq) = fpMB. (3.24) 

The constant fp is one of the key parameters of the heavy quark physics, just in 
the same way the constant is a key parameter of the soft pion physics. Below 
we will show that in the limit niq —>■ oo all three constants, fp, fv and f'p, coincide 

— 1/2 

with each other and scale as ruq modulo a weak logarithmic dependence on mq. 
(Needless to say, that both masses, Mp and My also coincide in this limit.) For 
definiteness let us consider f'p. Two other constants can be treated in a similar 
manner. The subscripts will be omitted in the remainder of this section to avoid 
overloaded expressions. 

Start from the two-point function 

A{k) = I J e^'=^rf^x(T{g(x)275g(x) m^^Qm) , (3.25) 
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where k is the external momentum. The two-point function A{k) develops a pole 
at k'^ = M^, the position of the ground state pseudoscalar, 



•^{k) = + excitations, (3.26) 

Of course, the currents produce from the vacuum not only the ground state mesons 
but also all excitations in the given channel. It is clear that to isolate the lowest- 
lying pole we should keep k"^ close to M^. Keeping in mind Eq. (|1.26| ) it is natural 
to represent k as 

k^ = {mQ + 6,0,0,0} 

where e scales like Aqcd while mq 00. With this parametrization of k^ we merely 
separate the mechanical (uninteresting) part of the momentum. The pole is achieved 
at e = A; near the pole 

The value of the coupling constant is obtained by amputating the pole. 

Let us now examine the theoretical expression for the same two-point function. 
In the background field technique (which is already pretty familiar, right?) we write 

A{k) = tTi ^^75— -1 I . (3.29) 

[ y -rriq l/+f-mQj 

Superficially this expression looks the same as if the quarks were treated as free; 
they are not, however; the coupling to the background field is reflected in the fact 
that Vfj, is the momentum operator, not just a c- number four- vector. 

Now we will take advantage of the fact that mg — > 00. As usual, we close our 
eyes on any possible hard contributions, assuming that V, the momentum operator 
of the light quark, is soft, i.e. does not scale with mg in the large mass limit but, 
rather V ~ Aqcd- (This is the reason, by the way, why the large external momentum 
k was directed through the heavy quark line in Eq. (|3.29|) .) Intuitively it is clear 



that the hard components of V should be irrelevant for the lowest-lying state whose 
"excitation energy" measured from itlq is of order Aqcd- 

If V is soft and tuq 00 the heavy quark Green function in the leading approx- 
imation takes the form 

1 ^l^+y + rnQ) ^ 



/+f-mQ ' {k + Vy -ml + {i/2)aG 

70+1 1 



:^ (3-30) 
/ n 



42 



where in the second hne all l/mq terms are omitted. No explicit mq dependence is 
left! This means that ^(e) scales as mg. Equation ( |3.28|) immediately implies then 
that / scales as 

/~mQ^/^ (3.31) 

Equation ( p.30|) for the heavy quark Green function in the limit mg ^ oo is in 
one-to-one correspondence with the leading term QTiQ{l + 7o)/2(5 in the low-energy 
Lagrangian ( |1.11| ). The analysis of the scaling law of the coupling constants pre- 
sented above is a simplified version of that carried out many years ago by Shuryak 
0. Later it was established that the power dependence on mq in Eq. ( |3.31| ) is sup- 
plemented by a logarithmic dependence appearing due to the hard gluon exchanges 



A few words about the numerical value of the coupling constants. Although 
in principle /d and fs are measurable experimentally, practically it is a very hard 
measurement, especially for B. No experimental number for exists so far. Its 
value was estimated in the QCD sum rules and on lattices more than once. Leaving 
aside a dramatic evolution of the issue I will say only that the recent and most 
reliable results cluster around 160 MeV both, in the sum rules [0, ^ and in the 



lattice calculations It is curious to note that in {rrihY^'^fB the preasymptotic 



l/mq correction turned out to be unexpectedly large and negative [^, ^ ^ ^ 



at the same time in {mbY^'^f^ the preasymptotic l/mq correction is much more 



modest 47, 48 



3.5 Proof of the Isgur-Wise formula 



I return to my promise to prove the Isgur-Wise formula ( p..25|) . Consider the three 



point function depicted on Fig. 6. The sides of the triangle are the Green functions 
of the quarks in the background gluon field. The reduction theorems tell us that in 
order to get the transition amplitudes {Hc\crb\Hf,) from this three-point function we 
must "amputate" it: multiply by (p^ — M|) and (p'^ — M|,), tending to M| andp'^ 
to M^). This singles out the meson states we want to pick up. For the vector mesons 
we must also multiply the three-point functions by its polarization vector e^. The 
last step necessary for amputation is dividing by the coupling constants (residues) 
connecting the currents bij^c and b'j^c to the respective mesons. If the currents are 
normalized appropriately - and we will always do that - the corresponding coupling 
constants in the pseudoscalar and vector channels are the same, fM (see Sect. 
3.4). It is convenient to combine the pseudoscalar and vector channels together by 
introducing the currents 

J = Bbi^^q + BfJ)'jf,q and J' = D*qi'j5C + D*^q-^^c (3.32) 

where B and D are external constants marking the annihilation and creation of the 
initial and final -B's and -D's {B^ and denote the polarization vectors of B* and 
D* , respectively). 
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The expression for the three-point function of Fig. 6 takes the form 



W'tTT^ ^W^] (3.33) 

I J/ + r -fric - nib y -rriq ) 



where the matrices A4' and A4 are introduced in Eq. ( |1.24 ). In Sect. 3.4 the 



niQ — * oo hmit of the quark propagator was obtained in the rest frame. Here 
we have two heavy flavor states, initial and final, and both can not be at rest 
simultaneously. Therefore, we need the very same propagator in the arbitrary frame. 
Let = rriQV^ + e^. Then a trivial generalization of Eq. ( |3.30| ) is 



j/+f-mQ ^ ' ij> + Vf -ml + (i/2)crG 

^+1 1 



(3.34) 



2 {e + V)v 

Using this propagator in Eq. ( p.34| ) we rewrite the three-point function of Fig. 



6 as follows 



V 2 2 J 



^ (e' + V)v' (e + V)v [f - m) J_ ' ^^'^^^ 

The expression in the braces is independent of the heavy quark masses; moreover, 
it is proportional to the three-point function in the theory with the scalar heavy 
quarks considered in Sect. 1.4. As was explained there, in this theory in the limit 
mg — i> oo only one form factor survives. 

One subtle point deserves discussing here. When I speak about the heavy flavor 
mesons I keep in mind particles built from the heavy quark and a light antiquark, 
which is not always in line with the accepted nomenclature. Say, they call B meson 
a particle with the b antiquark, not quark. Since this distinction plays no role in 
my lectures I will continue to ignore this linguistic nuance referring to the bq states 
as B mesons. All equations presented above assume that the b quark in the initial 
state annihilates to produce the c quark in the flnal state. Simultaneously a light 
antiquark in the initial meson is annihilated and the same light antiquark reappears 
in the flnal meson. 

Return now to the model with spinless heavy quarks. The heavy flavor hadrons 
we now deal with are spin-1/2 baryons. More exactly, we have antihaijon in the 
initial state and anizbaryon in the flnal state. This means that near the mass shell 
the expression in the braces in Eq. (|3.35| ) takes the form 



^^oMo J_^J_^ m: (3.36) 

/ pa 
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the minus sign between the unit term and the v term in the density matrices is due to 
the fact that we deal with the antibaryons. Amputating the legs and combining Eqs. 
{ ^^ ) and we get the Isgur-Wise formula ( |05| ) since A^(-^+l) = (^+l)A4, 

and so on. 



3.6 The Bjorken sum rule and all that 

The Isgur-Wise function C,{y) carries information about the structure of the light 
cloud. Needless to say that the heavy quark expansion per se does not help to 
calculate this function. One has to rely on methods applicable in the strong coupling 
regime which are outside the scope of my lectures (QCD sum rules, lattices, ...). Still, 
some interesting and important relations emerge. Here we will discuss a sum rule 
for the slope of the Isgur-Wise function and related topics. 

We are already familiar with the sum rule technology in the heavy quark theory. 
In Sect 2.1 we dwelled on a simplified problem: inclusive decays of a spinless heavy 
quark Q into a lighter spinless quark q and a fictitious spin-zero photon (j). The 
"photon" was assumed to be on mass shell, = 0. The predictions obtained referred 
to the moments of the "photon" energy. Now you are mature enough to face actual 
problems from real life. We will concentrate on the decays of a 6 containing hadron 
into a c containing hadron plus the lepton pair li>. The four-momentum of the lepton 
pair is a free parameter, in particular, ^ 0. We can and will choose the value of 
q to our advantage. 

Consider a transition Hh He induced by some particular current, say, axial- 
vector. At zero recoil ^ = 1. In the SV limit where the velocity of the recoiling 
hadron is small 



m = 1 - P\y - 1) + ... = 1 - p'y + - (3-37) 

where y = vv', v is the spatial velocity of in the Hf, rest frame, and the slope 
parameter was introduced in Ref. ||32[. It plays the same role as, say, the charge 



radius of pions. 

To get relations involving we start from consideration of the transition oper- 
ator similar to that of Eq. ( p.l3| ). The expectation value of the transition operator 
over the B meson state yields the hadronic amplitude whose imaginary part is pro- 
portional to the probability of the inclusive decay B XJiy with the fixed value of 
q, the momentum carried away by the lepton pair lu. (Here denotes an inclusive 
hadronic state containing one c quark.) A new element compared to the toy model 
of Sect. 2.1 is the heavy quark spin. Another distinction is the fact that, to achieve 
the SV limit, we do not need now to assume that rric is close to rrib. In the semilep- 
tonic decay B XJi' one can fine-tune the lepton pair momentum in such a way 
that q^ is close to its maximal value, g^g^^ — {^b ~ M^)'^; then the c containing 
hadronic state produced is almost at rest, and we are in the SV limit even though 
the charmed quark is significantly lighter than the b quark. In other words, for such 
values of g^ the c quark is always slow. 
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This transition operator describes the forward scattering of 5 to i? via inter- 
mediate states D* and "excitations" . (We focus for definiteness on the axial- vector 
current, cq^j^'-f^h. The vector current can be treated in a similar way.) The excitations 
can include, for instance, 0*7111, and so on. In general, all intermediate states except 
the lowest-lying D* will be referred to as excitations. The transition operator 



T 



I J ci^a;e-^^^T{6(x)7M75c(x) , cj,j,b} (3.38) 



in the Born approximation is given by the diagram of Fig. 1. The hadronic amplitude 
obtained by averaging T^j^ over the B meson state, 

huu = -^{B\f.jB) (3.39) 



contains various kinematical factors. In the general case the hadronic tensor h 



consists of five different covariant structures |3^, |33[| : 



-hg^^ + h2V^Vy - ihef^^a/sVaq/s + h^q^q^ + h^{q^v^ + q^v^). (3.40) 



Moreover, the invariant hadronic functions hi to /15 depend on two variables, go cind 
g^, or go and |g |. For g* = only one variable survives, and only two of five tensor 
structures in h^^ are independent. 

Each of these hadronic invariant functions satisfies a dispersion relation in go, 

= ^ / + polynomial (3.41) 

271 J go - go 

where Wi are observable structure functions, 

Wi = 2 Im hi . 

For our purposes it is quite sufficient to analyze only one function, namely, 
hi. Moreover, for the time being we will disregard all non-perturbative corrections 
^(^qcd) which means that operators in the expansion of T^i, other than hh can 
be neglected, and the B meson expectation value of hh can be replaced by unity. 
Calculating hi in this approximation is a trivial problem ( it was a part of the 
exercise suggested in Sect. 2.1). Specifically, 



h^^ = {mi, + me - go) - + O(Aqcd) (3-42) 



1 
— 

z 

where 

z = {m^-q^-E,){m,-qo + E,), Ej = ml + ^. (3.43) 

I remind that g is assumed to be fixed, and Aqcd -C |g| <^ Md, so that actually I 
will expand in g keeping only the terms up to second order. It is convenient to shift 
go by introducing a new variable 

e = qomax - go (3.44) 
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where 

qomax = Mb- Ed* , Ed* = Md* + • (3.45) 

ZIVI D* 

When e is real and positive we are on the physical cut where the actual intermediate 
states (e.g. D*) are produced. Here the imaginary part of hi is given by the "elastic" 
contribution of D* plus inelastic excitations. For negative e we are below the cut. 
The result for hi above can be trusted if — e 3> Aqcd since the expansion actually 
runs in Aqcd/c. The expansion in the inverse heavy quark mass also requires, of 
course, that |e| ^ Tric^i,. A bridge between the physical domain of positive e and 
the Euclidean domain of negative e where the calculation is done is provided by the 
dispersion relations. 

At the next stage the amplitude hi is expanded in powers of Aqcd/c and e/m^^c- 
Polynomials in e can be discarded since they have no imaginary part. We are 
interested only in negative powers of e. The coefficients in front of 1/e" are related, 
through dispersion relations, to the integrals over the imaginary part of hi with the 
weight functions proportional to the excitation energy to the power n — 1. Indeed, 

-hi{e, q ) = — de — = 

2tt J e — e 

Thus, our immediate task is to built the 1/e expansion from the amplitude ( p.42|) . 

The theoretical expression for the amplitude hi above knows nothing, of course, 
about the meson masses; it contains only the quark masses. Correspondingly, it is 
very convenient to build first the expansion of hi in an auxiliary quantity, 

e„ = nib - Ec- qo, Ec = mc + - — . (3.47) 

Then, if necessary, we reexpress the expansion obtained in this way in terms of e, 

1 = 1 + ^ + ... (3.48) 

The difference between and e is (9 (Aqcd ■ (f/^b,c) and O^AQQj^/mb^c)- 

After these introductory remarks, assembling all information in our disposal, we 

get 

-/.,= ^l--j- + A-- + ... (3.49) 

plus terms of higher order in (f or Aqcd- In deriving Eq. ( |3.49| ) I used the fact that 
e — eg = Aif /2 plus terms of higher order in or Aqcd- 

This completes the theoretical aspect of the calculation. The coefficients in front 
of 1/e and 1/e^ in hi are known; Eq. ( |3.46|) tells us that these coefficients are equal to 
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integrals over wi, the spectral density. So, what remains to be done is to express the 
spectral density in terms of the contribution coming from the physical intermediate 
states. Let us assume for simplicity that the spectrum of the intermediate states is 
discrete. Denote the mass of the z-th state by Mj and the energy by Ef, the lowest 
lying meson, D*, corresponds to i = 0. Then the propagator of the i-th meson 

1 



{Mb - qo - E,){Mb - qo + E,) 
at positive e has the imaginary part 

{2E,)-'7T6{e-6,) 

where 6i is the excitation energy (including the corresponding kinetic energy), 

Si = Ei — Ed* . 

For the "elastic" B — >■ D* transition 6q vanishes, of course. 

Now it is rather obvious that the structure function wi reduces to 



1/ 



|2 



«^i(^) = E^^2vr5(e-5.), (3.50) 

where the sum runs over all possible final hadronic states, the term with i = 
corresponds to the "elastic" transition B D* while i = 1,2, . . . represent excited 
states with the energies Ei = Mi + f/{2Mi); furthermore, l/s^ip looks like the 
square of a form factor. Strictly speaking l/^^jp is not exactly the square of a 
form factor; rather this is the (appropriately normalized) contribution to the given 
structure function coming from the multiplet of the degenerate states which includes 
summation over spin states as well. By appropriate normalization I mean that 
we routinely insert the normalization factor {2Mb)~^- In the particular example 
considered (the axial-vector current) D is not produced in the elastic transition, so 
that in the elastic part one needs to sum only over polarizations of D*. Say, at zero 
recoil fs^D* = ^/2Mj:,*Fb-^d* where Fb^d* is the B —>■ D* form factor at zero 
recoil, see Eq. ( [3.58| ) below. 

Let us examine in more detail the elastic contribution, i = 0. The form factor of 
the B D* transition generated by the axial- vector current is given in Eq. ( |1.21| ). 
Using this expression we readily obtain 

i2En*r\fB^nA' = ^ (^)\k)P - 1 - p'v' ■ (3-51) 

Comparing the 1/e coefficient in the dispersion representation ( p.46| ) with that 
of Eq. (|3.49| ) we conclude that 



+ (3.52, 
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which imphes, in turn 

1 °° I f |2 

In Sect. 1.3 we learnt that at zero recoil (i.e. v = 0) only the elastic transition sur- 
vives. As a consequence of the heavy quark symmetry for all non-elastic transitions 
iFg^jp ~ -iP. The ratio \fB^i\'^/v^ stays finite in the limit of small v. Eq. ( p.53|) 
is the Bjorken sum rule proper. Since the contribution of the excited states on the 
right-hand side is obviously positive it tells us, in particular, that > 1/4. This 
inequality is not very informative, though, since both, the QCD sum rule [^, |5^ 
and lattice calculations indicate that is only slightly less than unity, perhaps, 
close to 0.8. 

Leaving technicalities aside let me summarize the physical meaning of the result 
obtained. The coefficient in front of 1/e in Eq. ( |3.49|) does not contain Aqcd- This 
means that we calculate the probability of the decay b — c"Zi/" with given value of 
iP merely in the parton model; this probability is equal to that of the physical decay 
B Xc"/z/"; the latter is comprised of the elastic transition B 0*^^11/" and the 
transition of B into excitations. (The quotation marks are used to emphasize the 
fact that the decays that are measured are induced by both, the axial-vector and 
vector, currents, while we focus now only on the transitions induced by the axial- 
vector current.) All probabilities of production of the excited states are proportional 
to "iP (at small iP), and so is the part of the elastic transition containing p^. The 
sum of these two contributions must coincide with the "iP term obtained in the 
parton model. The very same analysis, by the way, presents a proof of the fact that 
C,{y = 1) = 1. (Do you see this?) 

Now, we make the next step, proceeding to the 1/e^ terms. The 1/e^ term in 
Eq. ( |3.49| ) is proportional to A, hence the result we are about to get evidently goes 
beyond the parton model. Combining Eq. ( p.49|) with the dispersion representation 
we find 



A- = V 5, = V (M, -Md,). (3.54) 

The sum runs not from zero to infinity but from 1 to infinity since = 0. Moreover, 
since all iFg^jp for z = 1, 2, ... are proportional to tP, and we are interested only in 
the iP term, it is legitimate to substitute, as I did, Ei by Mj and 6i by Mj — M^* . 
Eq. ( p. 541 ) is the optical (or Voloshin's) sum rule; superficially it looks the same 



as in the toy model of Sect. 2.1. Please, remember this sum rule - it gives a 
unique opportunity to measure A, one of the key parameters of the heavy quark 
theory. To this end one has to measure the inelastic transition probabilities in the 
semileptonic decays B XJiy in the SV limit. This is a difficult measurement, 
but not impossible, at least in principle. Before venturing into this noble task - 
extraction of A from experimental data - I must warn you that acceptable accuracy 
can be achieved only provided that the perturbative corrections (hard gluons) as 
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well as nonperturbative ones, of the next order in Aqcd, are included in the sum 
rules. We will briefly discuss the impact of the perturbative corrections in Lecture 
5. 

Those who are anxious to get something practical from the optical sum rule in 
the absence of the necessary measurements should not be discouraged. We can still 
get a lower bound on A. Indeed, let us rewrite Eq. ( |3.54| ) as follows: 



OO I £ |2 OO I J? |2 

Ay - g ^ (M. - M„.) + g ^ (M. - MO ^ 

(Afi - Mo.) (/ - \) if + E iyH^ {M. - Ml) ■ (3,55) 



where Eq. ( p.53|) is substituted. Since the second term on the right-hand side is 



obviously positive we conclude that 

A > 2(Mi - Md*) (p'^ ~\) ^ (3.56) 

where Mi is the mass of the first excited resonance with the quantum numbers of 
D* {Mi-Md* -0.5 GeV). 

Following the same line of reasoning one can derive the "third" sum rule relat- 
ing to an appropriately weighted sum over excitations ||5^. The corresponding 
inequality analogous to Eq. ( p.56| ) takes the form 



fil > 3(Mi - Md'Y (^P^ ~\) ^ • (3-^'^) 

3.7 Deviations of the B ^ form factor from unity at zero 
recoil 

The heavy quark theory began from the observation that the B ^ D* axial- vector 
form factor at zero recoil is exactly unity in the limit rrib oo , rric — > oo , mh/mc 
arbitrary, see Sect. 1.3. This is purely a symmetry statement, as usual, dynamics 
resides in the corrections. In this section we will discuss deviations from unity. 

When the heavy quark mass is infinite it is nailed at the origin, both in the initial 
B meson and in the final D*. The light cloud then does not notice the replacement of 
one quark by another, the overlap is unity. If we make the quark masses finite they 
start jiggling inside the mesons, and this motion is different in B and D* since the 
heavy quark velocities are different. On top of this the difference in the relative spin 
orientations of the heavy quarks and light clouds shows up. These two effects lead to 
deviations from unity. At a heuristic level there is no doubt that the deviations are 
of order of (i) the square of the characteristic heavy quark momentum {p itself can 
not enter since there is no preferred orientation) or (ii) chromomagnetic correlation 
aB. In both cases dimensional arguments prompt us that the deviation from unity 
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at zero recoil is proportional to l/m^^; linear effects in l/rric^i, are absent. The 
assertion was first formulated in Ref. [|15] and was cast in the form of a theorem 



(Luke's theorem) in Ref. |2T[]. The proof presented below is abstracted from the 
recent work [Q. 

Let us define the B D* form factor at zero recoil as follows 



(D*|c7^75&|5) = i^AMbM*^Fb^d*DI , (3.58) 

to be compared with Eq. (|L18| ). Conceptually our present derivation is very close 
to that leading to the Bjorken and Voloshin sum rules (Sect. 3.6). We will again 
consider the transition operator induced by the axial-vector current limiting our- 
selves to the spatial components of the current. Technically it is simultaneously 
simpler and more involved. Simpler - because at the point of zero recoil one must 
put g*= 0, so that kinematics is trivial. In particular, from the very beginning only 
one structure {hi) survives in the general decomposition ( p.40|) . The calculation is 
more complicated on the other hand since now one has to keep track of terms of 
order Aqcj-,. Those of order Aqcd are simply absent! 
The quantity e is defined now as 

e = Mb - Md, - go (3.59) 

and we continue to assume that Aqcd <^ |e| ^ m^fe and continue to examine our old 
acquaintance, hi, expanded in powers of 1/e and e/rricfi. The result of a relatively 
simple calculation (which the reader is encouraged to do) is 

- /ii = {1 - A} i + 0(A3jcd)^ + - (3.60) 

where 



A = + ^ ^ _ + _ + 

3 4 ml Smcmi,^ 



An explanatory remark is in order here concerning the term in hi. The 
theoretical expression for hi, as it naturally emerges from the computations, depends 
on tg, not on e where is the energy measured from the "quark" threshold, see Eq. 



( p. 47] ). In the kinematics at hand, when we are at the point of zero recoil, 

e- tg = Mb - Md* - {mi, - rric) = 



where I invoked Eq. ( 3.22|) . (Can you figure out why the coefficients in this ex- 



pression and in Eq. ( p.61| ) are different? Hint: The parameters and /i^ in Eq. 



( p.61| ) are defined as the expectation values of the corresponding operators over the 
pseudo scalar meson. This is not the case in Eq. (|3.22|) .) 
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We first expand —hi in l/e^ and then pass to the physical variable e and rearrange 
the expansion. In the l/e^ expansion the corrections of order O(Aqcd) are absent 
from the very beginning - an obvious fact hardly requiring further comments - while 
the term 0{Aqqy)/^1) does appear explicitly. This term, however, is killed in passing 
from 1/eq to 1/e, see Eq. ( ^.61| ). 



That is why I assert that the coefficient in front of is 0{Aqqy))- Moreover, 
as we will see shortly this coefficient in Eq. ( p.60|) is positive. In principle, it 



is calculable (more exactly, expressible in terms of several new phenomenological 
parameters) but this will be of no concern to us in this lecture. 

Repeating, step by step, the derivation of Sect. 3.6 we conclude that 



B^D* 



E \FB^e.cu\^ = 1 - A (3.62) 



=1,2,, 



and 

E \FB^e.cit\^ {Mi -Md^) = {A^cd)- (3-63) 

i=l,2,... 

The latter sum rule, by the way, is the reason why we know that the coefficient 
0{Aqqj)) in front of is positive - the left-hand side of the sum rule is obviously 
positive-definite. These two relations, taken together, plus positivity of iFg^jp, 
imply that \Fb^d*\'^ is limited from below and from above, 

1 - A - 1^^^%^ < \Fb^dA' < 1 - a (3.64) 

where Mi is the mass of the first excited state produced by the axial- vector current, 
Mi-Md* =0(Aqcd). 

Not only is it seen that the deviation of |Fb^d*P from unity starts from terms 
scaling like with no l/rric^b corrections, but we understand now the reasons 

lying behind this remarkable fact. Moreover, we have an idea of how large the 
actual deviations are since Eq. (|3.64D establishes a lower limit for these deviations 



in terms of the parameter A which is determined numerically rather well. In this 
aspect the derivation I present here goes beyond the more conventional analysis of 



Ref. 1^3], 0. The reader is nevertheless advised to consult the latter works to get 
a broader perspective of the heavy quark theory - the more approaches you master 
the better for you. 

Qualitatively it is quite clear why the deviation of {Fb^d-A"^ from unity starts 
from Aqcj-j/m^^. Indeed, let us return to the Bjorken formula ( |3.37| ). In this formula 
it is assumed v ^ Aqcd/^^q so that actually we do not distinguish between the 
velocity of the recoiling final heavy hadron and that of the final quark. At zero 
recoil the heavy hadron is nailed, but not the heavy quark. The latter experiences 
a primordial motion inside the nailed hadron, with the velocity -iP ~ AqQ^/mQ. So, 
a reasonable guess would be to extrapolate Eq. (|3.37|) down to v"^ ~ ^qcd/''^q- ^s 
we see, this guess works. 
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It is worth emphasizing that our analysis need not be confined to the transitions 
induced by the spatial components of the axial- vector current. We could consider 
the temporal components, or vector currents, or something else. Each time we get 
additional information. For instance, from the transition operator induced by the 
vector currents we get a sum rule proving the inequality //^ > /Xq obtained in Sect. 
3.2 from a quantum-mechanical argument. 
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4 Lecture 4. Theory of the Line Shape 



In this lecture I will discuss one of the most interesting and practically important 
applications of the heavy quark theory, the spectra in the end point domain in the 
inclusive decays. Inclusive weak decays of heavy flavors, in particular, semileptonic 
decays, are close relatives of famous deep inelastic scattering - the processes where 
a highly virtual photon scatters off nucleons to produce an inclusive multiparticle 
state. The latter are related to the former via channel crossing. Deep inelastic 
lepton-nucleon scattering was in the focus of theoretical activity in the late sixties 
and the beginning of seventies and was instrumental in discovering and developing 
QCD [Q. It is thus quite surprising that for a long time there were hardly any 
attempts to treat the beauty decays in QCD proper along essentially the same lines 
as it was done in deep inelastic scattering. Realization of the idea that the l/mq 
expansion in the theory of the line shape can play the same role as the twist 1/Q^ 
expansion in DIS came with the 20 years delay ^ - I see absolutely no 

reasons why the corresponding theory was worked out only recently and not 20 
years ago. 

The theory of the line shape in QCD resembles that of the Mossbauer effect. 
To explain what I mean it is convenient to consider, for definiteness, the transition 
B — »• where Xg denotes the inclusive hadronic state with the s quark. This 
decay has been recently observed experimentally. (Description oi B ^ Xqlu is 
conceptually similar but is more technically involved). 

Again, to avoid inessential technicalities I will neglect the quark and photon 
spins. So we will consider the transition Q — g0 where all fields Q, q and are 
spinless. Thus, to begin with, we will limit ourselves to the toy model described in 
Sect. 2.1, see Eq. (O). The mass of the final quark niq will be treated as a free 
parameter which can vary from zero almost up to mg. For our approach to be valid 
we still need that Am = niq — ruq ^ Aqcd although the mass difference may be 
small compared to the quark masses. 

To warm up we will put the final quark mass to zero. At the level of the free 
quark decay the photon energy is then fixed by the two-body kinematics of the decay 
Q ~^ Q'Pi namely, = mQ/2. In other words, in the rest frame of the decaying 
Q quark the photon energy spectrum is a monochromatic line at = mQ/2 (Fig. 
7). On the other hand, in the actual hadronic decays Hq — >■ Xq(j) the kinematical 
boundary of the spectrum lies at Mhq/2. Moreover, due to multiparticle final states 
(which are, of course, present at the level of the hadronic decays) the "photon" line 
will be smeared. In particular, the window - a gap between mQ/2 and Mhq/2 - will 
be closed (Fig. 7). There are two mechanisms smearing the monochromatic line of 
the free-quark decay. The first is purely perturbative: the final quark q can shake off 
a hard gluon, thus leading to the three-body kinematics. This mechanism tends to 
diminish the photon energy and may be important at £^ < mQ/2. We will defer its 
discussion till later times. The second mechanism is due to the "primordial" motion 
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of the heavy quark Q inside Hq and is non-perturbative. Even if the decaying B 
is nailed at the origin so that its velocity vanishes, the b quark moves inside the 
hght cloud, its momentum being of order Aqcd- This is the QCD analog of the 
Fermi motion of the nucleons in the nuclei. It is quite clear that this motion affects 
the decay spectra. Say, if the "primordial" heavy quark momentum is parallel to 
that of the photon, the photon produced gets more energy, and vice versa, for the 
antiparallel momenta it gets less. It is quite clear that this effect is preasymptotic 
(suppressed by inverse powers of nib): while typical energies of the decay products 
are of order rrih a shift due to the heavy quark motion is of order Aqcd- 

Only the second mechanism will be of interest for us in this lecture. 

The window (i.e. the domain kinematically inaccessible for free quarks) plus the 
adjacent domain below the window, of width several units xAqcd, taken together, 
form what is called the end point domain. Below I will outline the main elements of 
the theory allowing one to translate an intuitive picture of the Q quark primordial 
motion inside Hq in QCD-based predictions for the spectrum in the end point 
domain . The spectrum below the end point domain is the realm of the perturbative 
physics (hard gluon emission). 



4.1 Formalism 

Let us return back to Sect. 2.1 and consider the transition operator defined there. 
Since we are interested in the energy spectrum the "photon" momentum q must 
be fixed. Let us assemble Eqs. ( |2.14| ) and ( |2.23| ) together. For convenience I will 



reproduce the result here again, taking into account the fact that now is assumed 
to vanish. 



dE Mhq V 2 / AruQ \ 2 / 24 V 2 / J 

(4.1) 

If in the leading approximation the spectrum is just a delta function, the corrections 
are more and more singular! The higher the correction the stronger the singularity. 
Nonsense? No, this was to be expected: the width of the line in the transition 
Hq — > Xq(f) is of order Aqcd- We expand in the powers of Aqcd/'^q; hence we must 
expect the enhancement of the singularities in each successive order. Equation ( [4.1|) 
gives all terms up to Aqcd- It is clear that to describe the shape of the line one 
needs to sum up the infinite number of terms in this expansion. 

Then in the approximation of Fig. 1 (no hard gluon exchanges) the transition 
operator is given by Eq. ( p.l5|) with niq set equal to zero. To construct the operator 



product expansion to all orders we observe that the momentum operator n corre- 
sponding to the residual motion of the heavy quark is ~ Aqcd and the expansion 
in ir/k is possible. Unlike the problem of the total widths, however, in the end 

^Basics of the theory of the line shape were worked out in Refs. j5^, Further crucial 

steps were undertaken in Refs. |60, 6^, In my presentation I follow mainly Bigi et al. [ p9[ . 



55 



point domain k"^ is anomalously small, the expansion parameter is of order unity, 
and there exists an infinite set of terms where all terms are of the same order of 
magnitude. 

To elucidate this statement let us examine different terms in the denominator of 
the propagator, 

k"^ + 2'nk + TT^ 

In the end point domain 

|E- (l/2)mQ| ~ A (4.2) 
It is quite trivial to find that in this domain 

fco ~ \k\ ~ mQ/2, fc^ ~ rnqK] 

in particular, at the kinematical boundary (for the maximal value of the "photon" 
energy) k^ = Mq/2 and /c^ = —itlqA. Hence, in the end point domain 

k'^ r-u kn ^ TT^. 

In other words, when one expands the propagator of the final quark in the transition 
operator, 

in TT, in the leading approximation all terms (2A;7r/fc^)" must be taken into account 
while terms containing vr^ can be omitted. The first subleading correction would 
contain one tt^ and arbitrary number of 2/c7r's, etc. 

Thus, in this problem it is twist of the operators (= dimension - Lorentz spin) in 
the operator product expansion, not their dimension, that counts. For connoisseurs 
I will add that this aspect makes the theory of the line shape in the end point 
domain akin to that of deep inelastic scattering (DIS). Keeping only those terms 
in the expansion that do not vanish in the limit rriQ oo (analogs of the twist-2 
operators in DIS) we get the following series for the transition operator 

-| oo y 9 \ " 

2' = -pE(-;^) A^'^^..A^'^"(g7r^....vr^„Q- traces). (4.4) 

Traces are subtracted by hand since they are irrelevant anyway; their contribution 
is suppressed as k'^Ti'^ / {knY ~ Aqcd/^^q to a positive power. Another way to make 
the same statement is to say that in Eq. ( |4.4| ) the four-vector k can be considered 
as light-like, /c^ = 0. 



4.2 The light cone distribution function 

After the transition operator is built the next step is averaging of T over the hadronic 
state Hq. Using only the general arguments of the Lorentz covariance one can write 

{Hq\Qtx^^...'k^^^Q - tracesliJg) = a„A"(f^^...t;^„ - traces) (4.5) 
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where a„ are constants parametrizing the matrix elements. Their physical meaning 
will become clear momentarily. Right now it is worth noting that the term with 
n = 1 drops out (oi = 0). Indeed, {Hq\Qt:Q\Hq) is obviously zero for spinless Hq 
while ttq through the equation of motion reduces to 7f^/(2mQ) and is of the next 
order in l/mq. Disappearance of Qtc^Q means that there is a gap in dimensions of 
the relevant operators. 

Let us write a„'s as moments of some function F{x), 

an = J dxxTix). (4.6) 

Then, F{x) is nothing else than the primordial line-shape function! (That is to 
say, F{x) determines the shape of the line before it is deformed by hard gluon 
radiation; this latter deformation is controllable by perturbative QCD). The variable 
X is related to the photon energy. 




If this interpretation is accepted - and I will prove that it is correct - it immediately 
implies that (i) F{x) > 0, (ii) the upper limit of integration in Eq. ( [4.6[ ) is 1, (iii) 
F{x) exponentially falls off at negative values of x so that practically the integration 
domain in Eq. ( [4.6| ) is limited from below at —xq where Xq is a positive number of 
order unity. 

To see that the above statement is indeed valid we substitute Eqs. ( |4.5|) , ( |4.6|) 
in T, 

(W.}--^i:l''ymr{-^J. (4.7) 

and sum up the series. The iO regularization will prompt us how to take the imagi- 
nary part at the very end. In this way we arrive at 

^ = /^^'^'FTi^krio = (2M)ronx). (4.8) 

where the variable x = k'^/{2Akv) was written out above in terms of the "photon" 
energy, and Fq is the total decay width in the parton approximation. Corrections 
to Eq. ( [4. 81) are of order A/mg. 

Thus, we succeeded in getting the desired smearing: the monochromatic line of 
the parton approximation is replaced by a finite size line whose width is of order A. 
The pre- asymptotic effect we deal with is linear in Aqcd/'^q- 

At this point you might ask me how this could possibly happen. We have already 
learnt that there is a gap in dimensions of the operators in the expansion - no 
operators of dimension 4 exist - and the correction to QQ is also quadratic in 1/ mq 
(the CGG/BUV theorem). There are no miracles - the occurrence of the effect 
linear in Aqcd/'^q became possible due to the summation of the infinite series in 
Eq. (|4.4|) ; no individual term in this series gives rise to Aqcb/^q- 
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To avoid misunderstanding it is worth explicitly stating that the primordial 
distribution function F{x) is not calculated; rather F{x) is related to the light-cone 
distribution function of the heavy quark inside Hq, namely {Q{rni)'^Q) , = 0, or 
more explicitly 

F{x) oc j dte'^'^{HQ\Q{x = 0)e-^/o "^("^)'='^Q(a;^ = (4-9) 
where n is a light-like vector^ 

n^ = (1,0,0,1). 

Unfortunately, this primordial function is not the one that will be eventually mea- 
sured from dr/dE; the actual measured line shape will be essentially deformed by 
radiation of hard gluons. I will say a few words about this in Lecture 5. 

The primordial distribution function F{x) which we defined here can be called 
the light- cone distribution function. This is clear from the expression ( |4.9| ) which 
has a very transparent physical meaning. The quark q produced is massless and, 
therefore, propagates along the light cone from the point of emission to the point of 
absorption in the transition operator defining the distribution function. 

If we looked at the physical line shape sketched on Fig. 7 more attentively, 
through a microscope, we would notice that a smooth curve is obtained as a result 
of adding up many channels, specific decay modes. A typical interval in E that 
contains already enough channels to yield a smooth curve after summation is ~ 
^qcb/^Q- Roughly one can say that the spectrum of Fig. 7 covers altogether 
mg/AQCD resonance states produced in the Hq decays and composed of q plus the 
spectator (I keep in mind here that the final hadronic state is produced through 
decays of highly excited resonances, as in the multicolor QCD). These states span 
the window between mQ/2 and and the adjacent domain to the left of the 

maximum at E = mQ/2. 



4.3 Varying the mass of the final quark 

So far I was discussing the transition into a massless final quark. It is very interesting 
to trace what happens with the line shape and the primordial distribution function 
as the final quark mass rriq increases. 

Inspection of the transition operator shows that as long as <^ wqA nothing 
changes in our formulae at all in the leading-twist approximation. Since the char- 
acteristic values of k"^ in the end point domain are of order ttlqA and <^ rnQA 
one can merely neglect the final quark mass altogether. 

A more interesting regime is niq ~ (mgA)^/^. In this regime one can not neglect 



mz in the denominator. It is not difficult to see, however Il62|l, that, as in the massless 



10 



Similar light-cone distributions for light quarks are well known [M in the theory of deep 



inelastic scattering, see also |65 . 
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case, all traces can be neglected since k'^ir'^/lkir)'^ ~ Aqcd/'^q- This means that 
the same light-cone distribution function F{x) that emerged in the massless case 
describes the line shape if rUg ~ mgA as well. The only change that occurs is a shift 
of the end point spectrum, as a whole, to the left. Indeed, if previously the variable 
X was defined as {2/A)[E — (l/2)mQ], now when rriq 7^ 

x=j{E-Eo) 

where Eq = {2mQ)^^{mQ — m^). The maximum of the distribution, in particular, 
shifts from mQ/2 to Eq = mg/2 - C(A). 

What happens if one continues to increase m^? Increasing the quark mass further 
results in more drastic changes. The trace terms can not be omitted any more, and 
the light-cone function gives place to other distribution functions. This is obvious 
already from a simple kinematical argument. Indeed, with increasing the window 
shrinks. When we eventually come to the SV limit 

Aqcd < Am = mg - mq < mg^q 

it shrinks to zero. In this limit the photon energy in the two-body quark decay, 
Am(l + Am(2mQ)~^) differs from the maximal photon energy in the hadronic decay, 
AM(1 -|- AM{2Mq)~^), only by a tiny amount inversely proportional to mq (Am 
and AM stand for the quark and meson mass differences, respectively). 

Thus, the kinematical consideration prompts us that the line shape must es- 
sentially change. Anticipating the results of the calculation let me describe the 
situation pictorially. Simultaneously with the shrinkage of the window the peak 
becomes more asymmetric and develops a two-component structure (Fig. 8). The 
dominant component of the peak, on its right-hand side, becomes narrower and 
eventually collapses into a delta function when rriq becomes a finite fraction of rriQ. 
A shoulder develops on the left-hand side; the number of the hadronic states popu- 
lating the end point domain becomes smaller - instead of mg/AqcD states at ruq = 
we are speaking of just several states at rUq = a finite fraction of tjiq. When we 
approach the SV limit the height of the shoulder corresponding to the production of 
the excited states becomes very small, proportional to -y^ <^ 1 (Fig. 9). This is the 
end of the evolution - starting from the light-cone distribution function at m^ = 
we continuously pass to the temporal distribution function in the SV limit. It is the 
temporal distribution function that shapes the inelastic shoulder on Fig. 9. 

This rather sophisticated picture, hardly reproducible in naive quark models, 
emerges from the operator product expansion (in the leading approximation) if one 
follows along the same lines as previously. The transition operator T for rUq 7^ is 
given in Eq. (|2.20| ); I reproduce it here again for convenience. 
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Notice that 2mQno + tt^ acting on Q yields zero (the equation of motion) and in the 
SV hmit q must be treated small parameter, 



qo/mQ = E/niQ = < 1; 

V is the spatial velocity of the heavy quark produced. Although v is small the 
inclusive description is still valid provided that Am ^ Aqcd- 

In the zeroth order in q the only term surviving in the sum ( [4.10| ) is that with 
n = 0, and we are left with the single pole, the elastic contribution depicted on Fig. 
9. This is the extreme realization of the quark-hadron duality. The inclusive width is 
fully saturated by a single elastic peak. We have already discussed this phenomenon 
in Lecture 2. What might seem to be a miracle at first sight has a symmetry 
explanation - the phenomenon is explained by the heavy quark symmetry. The fact 
that the parton-model monochromatic line is a survivor of hadronization is akin to 
the Mossbauer effect. 

If terms (9(f^) are switched on the transition operator acquires an additional 
part, 

= 5«"°K^ I (4^)" Q--<--Q- (4.11) 

From this expression it is obvious that the shape of the v"^ shoulder is given by 
the temporal distribution function G{x) whose moments are introduced through the 
matrix elements 

{Hq\Q7t,7t^7t,Q\Hq) = A"+2 J rfxx"G(x). (4.12) 
Alternatively, G{x) can be written as a Wilson line along the time direction, 

G{x) oc J dte''''~^{HQ\Q{t = 0, f = 0)7r,e-'/o ^''W'^V,g(t, f = 0)|i/Q). (4.13) 

Intuitively it is quite clear why the light-cone distribution function gives place 
to the temporal one in the SV limit. Indeed, if the massless final quark propagates 
along the light-cone, for Am -C m the quark q is at rest in the rest frame of Q, i.e. 
propagates only in time. 

In terms of G{x) our prediction for the line shape following from Eq. ( [4.11|) takes 
the form 

dT 



dE"" 



1-^ 
3 



S{x)+ 




X 



G(x), (4.14) 



where x = {E — E„iax)/^- The t>^ corrections affect both, the elastic peak (they 
reduce the height of the peak) and the shoulder (they create the shoulder). The total 
decay rate stays intact, however: the suppression of the elastic peak is compensated 
by the integral over the inelastic contributions in the shoulder. This is the Bjorken 
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sum rule thoroughly considered in Lecture 3. It is important that we do not have to 
guess or make ad hoc assumptions - a situation typical for model-building - QCD 
itself tells us what distribution function enters in this or that case and in what 
particular way. 



4.4 Real QCD: Inclusive semileptonic decays 

From the analysis presented above the following remarkable fact should be clear. 
The very same primordial distribution functions that determine the line shape in 
the radiative transitions appear in the problem of the spectra in the semileptonic 
decays. In particular, in 6 — ulu we deal with F{x). 

Of course, kinematical conditions are different. Now the hadronic part of the 
process, B XJu inclusive decay, depends on two variables, for instance, go and 
g^, or go and |g|. The probability of the decay, in the free quark approximation, is 
proportional to S{mQ — go — |g|) [0- In other words, in this approximation only a 
line on the go, |g| plane is populated (Fig. 10). (I assume that we are not interested 
in the individual momenta of I and u and measure just the total momentum of the 
lepton pair. This is quite a fantastic formulation of the problem since experimentally 
the neutrino energy and momentum are not measured, of course; only the electron 
energy is usually measured. Nevermind, let us keep in mind a ^^edan/sen experiment.) 
The end point domain is defined now as a band whose width is several units xAqcd 
adjacent to the above quark line (Fig. 10). Needless to say that in the physical 
decay the whole large triangle is populated; the inner part of the triangle, to the 
left of the end point band, is due to the hard gluon emission. The smearing of the 
delta-like spectrum in the band is due to the primordial motion of b inside B, and 
is described by the light cone distribution. 

A trivial modification compared to Sect. 4.3 is the occurrence of several structure 
functions. All five structure functions are expressible, however, in terms of the 
same light cone primordial distribution function F{x) where, as previously, x = 
—A'^k"^ /2ko. Since go and g are independent variables in the case at hand 

X = -A-'ky2ko = -A-'k^iko + \k\) = -A-\mQ - go - |gl) (4.15) 

where in the denominator the difference between k^ and \k\ is neglected which is 
perfectly legitimate in the end point band. In this band k^ — \k\ = O(Aqcd), and the 
difference between fco and \k\ becomes important only at the level of the subleading 
twists which are not included anyway. 

Thus, we observe a scaling behavior: the structure functions that generally speak- 
ing could depend on two variables, go and |g|, actually depend only on the single 
light-cone combination ( |4.15| ). This is the analog of the Bjorken scaling in deep 



inelastic scattering! In the rest of the phase space, outside the end point band, the 
approximate equality k^ ^ \k\ is not valid, of course, and the above scaling is not 
going to take place. The primordial distribution falls off - presumably exponentially 
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- outside the end point band. The hard gluon emissions will populate the phase 
space outside this domain creating long logarithmic tails. The primordial part is 
buried under these tails. Therefore, outside the band one can not expect that the 
structure functions depend on the single combination go + |^ anyway. 

Guesses about a scaling behavior in the inclusive semileptonic decays are known 



in the literature |6§]. Now we are finally able to say for sure what sort of scaling 
takes place, where it is expected to hold and where and how it will be violated. 

I will not go into further details which are certainly important if one addresses 
the problem of extraction of V^b from experimental data. Some of them are discussed 
in the literature , others still have to be worked out. Applications of the theory to 
data analysis is a separate topic going beyond the scope of this lecture. 

What can be said about the hght cone distribution function This function 

depends on the structure of the light cloud of the B meson and, thus, belongs to the 
realm of the soft physics. The moments of this function are related to the expectation 
values of the operators Qtx^^...t[^^Q (see Eq. ( [4.5| )); in real QCD the properly 
normalized matrix elements on the left-hand side include the factor {2Mb)~^)- The 
knowledge of the infinite set of these expectation values would be equivalent to the 
knowledge of the structure of the light cloud. Needless to say that this is beyond our 
abilities at present. Still, we know a few first moments of F{x) and have a general 
idea of the shape of this function. It must be positive everywhere in the physical 
domain, vanish at x = 1 and have exponential fall-off at large negative x. The latter 
property ensures the existence of all moments. Moreover, 



a, = jdxF[x) = l, 
tti = J dxxF(x) = - 



and 



a, = j dxx^Fix) 



3A2 ■ 

Estimates of the third moment also exist in the literature ||59|, I can not dwell 
on this issue now and will only mention that 03 is constrained by exact inequalities. 



i.e. |6l 



1/1 If ly 

a2 < - + ]l--a3 , a3<--[a2--j . (4.16) 

To derive these inequalities one merely observes that for any t the integral from 
—00 to 1 over X over the function (1 — x){x — tYF{x) is positive; on the other hand 
this integral is a second order polynomial in t and, hence, its discriminant must be 
negative. 

A sketch of a function satisfying all these requirements is given on Fig. 11. 

A natural desire to extend the formalism described above to the semileptonic 
inclusive transitions h civ encounters serious technical difficulties. The essence 
of the problem is as follows. The final quark c can be treated as heavy, although 
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at the same time , -C m^. The ratio ml/ml ^ 0.07 is a small parameter 
while ml/{Amb) ~ 1. Under the circumstances the type of the distribution function 
describing the primordial motion of b inside B and determining the measurable 
structure functions Wi to will depend on the value of and the scaling property 
- dependence on one particular combination of variables - is lost. 

The Qq, \q\ plane is shown on Fig. 12. In the free quark approximation the 
transition probability is proportional to 5{mh — % — E^) where Eq = (m^ + (fY^'^, 
and all events are concentrated along the line indicated on Fig. 12. At the hadronic 
level the phase space consists of the full triangle, with one side curved. The end 
point band is also curved. 

The fact that one side of the triangle is distorted compared to 6 — >■ ulu is not 
crucial. What is important is the change of dynamics as we move from the upper 
left corner to the lower right one. In the case of 6 — * ulu moving along the end 
point band in this direction does not affect the measured structure functions (apart 
from the extreme domain of soft u - the exclusive resonance domain - where our 
description fails altogether). The situation is different in the b — *• diy transition. 

If ^ ml one recovers the same light-cone function F{x) as in the tran- 
sition b ulv OT b —>■ s'-f. Modifications are marginal. First, some extra terms 
explicitly proportional to mc/ m^ are generated in the structure functions due to the 
fact that y + mc replaces y in the numerator of the quark Green function. Moreover, 
if in the b ulv transitions the scaling variable in the end-point domain is 

x = A-^(go + |<?1-mb), (4.17) 

in the b civ transition it is shifted by a constant term of order 1, 

2 

— mi 
x = k~\qo+\q\-mb) + ^. (4.18) 

Am/, 



To see how this shift occurs and to reveal limitations of the approximation let 
us start from the parton model variable m^ — qo — E^.. In the limit |^^/m^ oo 
inside the end point band formally one may substitute E^ by 

2 2 

mt , ^ m: 



If ml = 0{Amb) formally one may discard the 0{ml/ml) correction. In this way 
we arrive at the scaling variable ( [4.181) . 



It is worth emphasizing that the occurrence of the light-cone distribution function 
in this regime, the same as in the 6 — > m transition, is a remarkable fact. Indeed, 
if we could examine the measured structure functions "in the microscope" in these 
two cases we would see that their microstructure is quite different. As was already 
mentioned, in the b u transition the end point band is saturated by the production 
of ~ mb/A states, with the spacing between the individual states of order A^/mfe. 
In the b —>■ c transition, even if we are in the upper left corner of the phase space 
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where the hght-cone distribution is relevant, the number of the states produced is 
~ mfy/rric ~ rric/A and the spacing is of order A^/mc. We deal with a much coarser 
structure in the latter case, and still all resonance contributions being summed up 
must add up to produce the light-cone distribution, formally the same one that is 
created by a much larger number of resonances in the b —>■ u transition. (Purely 
theoretically we can not predict fine grain versus coarse grain composition of the 
structure functions if we limit ourselves to the leading twist. Only analysis of all 
twists could resolve these details, would this analysis be possible). 

The word "formally" is used above three times, not accidentally. Practically in 
the b civ transition \q\ can never be much larger than m^. Indeed, the maximal 
value of corresponding to = is (m^ — m^)/2mb ~ 2 GeV, that is only 
~ 1.5mc. Therefore, only by stretching a point and only in a narrow domain near 

= 0, one can expect that the light cone function of the variable ( [4.18|) is, perhaps, 
more or less relevant. 

As |g| decreases and becomes less than rric (this regime takes place in a large 
part of the phase space) the light-cone distribution function becomes irrelevant. 
The measurable structure functions are determined by a different distribution - the 
light-like vector in Eq. (|4.9|) is replaced by 




and it becomes clear that the would be scaling variable x = A^^(go + Ec — mb) 
fails to represent all dependence of the structure functions on go aiid |g|. When 
approaches its maximal value, 

qL. = (Mb - Mof 



|g| tends to zero and we eventually approach the SV regime which I have already 
discussed, with fascination, in the toy example above. In the SV limit the velocity 
of He produced is small, and the structure functions probe the primordial motion 
described by the temporal distribution function G{x) where now 

X ^ A"^(go — Am), 

Am is the quark mass difference coinciding, to the leading order with Mb — Mb,. 

Thus, changing q^ from zero to q^^x results in an evolution of the distribution 
function appearing in theoretical formulae for dT{B XJu), from light-cone to 
temporal, through a series of intermediate distributions. The physical reason for this 
evolution is quite clear - what distribution function is actually measured depends 
on the parton-model velocity of the quark produced in the b decay. In the limiting 
cases of very large recoil and very small recoil the problem is solved in the sense 
that the structure functions are expressed in terms of the light-cone and temporal 
distribution functions, respectively. The intermediate case |g| ~ rric is not worked 
out in detail so far. It is beyond any doubt, however, that the parton-model type 
scaling will not take place. 
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5 Lecture 5. Including Hard Gluons. Generali- 
ties of the Operator Product Expansion. 

Finally the time comes when I can not ignore any more the existence of hard gluons. 
Hard gluons are mere nuisance from the point of view of the theory of hadrons since 
they play no, or very little, role in the structure of the low-lying hadronic states. Yet, 
if we want to go beyond purely academic exercises, however beautiful they might 
look, and descend down into a messy world of real hadronic physics, hard gluons 
can not be forgotten about since they "contaminate" nearly every experimentally 
measurable quantity. To make contact with the real world we have to consider 
interplay between the soft and hard physics. 

The hard gluons manifest themselves in many ways. They contribute to the 
coefficient functions in the effective Lagrangian (|1.3|) obtained by integrating out all 
degrees of freedom with the characteristic frequencies down to fi. They show up in 
the calculations of the total decay rates and spectra discussed in Lectures 3 and 4 
resulting in perturbative corrections which, in some instances, change the answer 
quite drastically. They result in the fact that all basic parameters of the heavy quark 
physics - the heavy quark mass. A, fi^. and so on - generally speaking, become fi 
dependent and can not be treated as universal constants. Here we will address some 
of these issues in brief. 



5.1 Calculation of the effective Lagrangian 

I have already started discussing this topic in Sect. 1.2. The original QCD La- 
grangian ( p. . 1|) is formulated at very short distances. In principle, it codes all in- 
formation necessary for calculation of all observable amplitudes. We just do the 
functional integral and ... Alas, there are very few functional integrals that can be 
calculated analytically; numerical evaluation on lattices may take years, and I even 
dare to assert that some amplitudes will never be calculated that way. So, we take 
the original Lagrangian and start evolving it down, integrating out all fluctuations 
with the frequencies fi < u < Mq where Mq is the original normalization point, and 
fi will be treated, for the time being, as a current parameter. In this way we get the 
Lagrangian which has the form 

C = J2Cn{Mo;fi)OM. (5.1) 

n 

The coefficient functions C„ represent the contribution of virtual momenta from fi to 
Mq. The operators On enjoy full rights of the Heisenberg operators with respect to 
all field fluctuations with frequencies less than fi. The sum in Eq. ( |5.1| ) is inflnite - it 
runs over all possible Lorentz singlet gauge invariant operators with the appropriate 
quantum numbers; for instance, if CP is conserved, only CP-even operators will 



appear in (|5.1|). If, say, the electromagnetic processes are included, the operators 



in the Lagrangian (|5.1|) may contain the photon and electron flelds, and so on. All 
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operators can be ordered according to their dimension; moreover, we can use the 
equations of motion stemming from the original QCD Lagrangian to get rid of some 
of the operators in the sum. Those operators that are reducible to full derivatives 
give vanishing contributions to the physical (on mass shell) matrix elements and can 
thus be discarded as well. 

If one just abstractly writes the expression one is free to take any value of 
fi] in particular, /i = would mean that everything is calculated and we have the full 
S matrix, all conceivable amplitudes, at our disposal. Nothing is left to be done. In 
this case Eq. ( |5.1|) is just a sum of all possible amplitudes. This sum then must be 
written in terms of the physical hadronic states, of course, not in terms of the quark 
and gluon operators since the latter degrees of freedom are simply non-existent at 
large distances. 

This is day-dreaming, of course. Needless to say that in our explicit calculation 
of the coefficient functions we have to stop somewhere, at such virtualities that the 
quark and gluon degrees of freedom are still relevant, and the coefficient functions 
Cn(Mo, fi) are still explicitly calculable. On the other hand, for obvious reasons it is 
highly desirable to have fi as low as possible. In the heavy quark theory there is an 
additional requirement that fi must be much less than mq. The process of calculat- 
ing the coefficients Cn{MQ,fi) is called matching in the more standard presentation 
of HQET. Actually we see that this procedure is nothing else than a generaliza- 
tion of Wilson's idea of the renormalization group and the (Wilsonean) operator 
product expansion. Using the standard OPE language has an evident advantage: 
all well-studied elements of the latter approach can be immediately adapted in the 
environment of the heavy quark expansions. In particular all parameters one can 
read off from the Lagrangian (|5.1| ) depend on fi (including, say, the heavy quark 
mass). Let us assume that is large enough so that Q;s(/i)/7r ^ 1, on the one 
hand, and small enough so that there is no large gap between Aqcd and fi. The 
possibility to make such a choice of n could not be anticipated apriori and is an 
extremely fortunate feature of QCD, a gift from the Gods. Quarks and gluons with 
the offshellness larger than /i chosen that way are called hard. 

Needless to say that the parameter n is in our minds, not in Nature. All ob- 
servable amplitudes must be fi independent. The /i dependence of the coefficient 
functions C„ must conspire with that of the matrix elements of the operators On in 
such a way as to ensure this /i independence of the physical amplitudes. 

What can be said about the calculation of the coefficients C„ ? Since fi is suffi- 
ciently large, see above, the main contribution comes from perturbation theory. We 
just draw all relevant Feynman graphs and calculate them, generating an expansion 
in as(/i) which for brevity I will denote by a^, with the argument omitted, 

Cn = ^aia[ . 
I 

Sometimes some graphs will contain not only powers of as(/i) but powers of ln(mQ//i). 
This happens if the anomalous dimension of the operator On is nonvanishing - quite 
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a typical situation - or if a part of a contribution to C„ comes from characteristic 
momenta of order rriQ and is, thus, expressible in terms of ^^(mQ), and we rewrite 
it in terms of ^^(/i). Nevermind, this is a trivial technicality. You are supposed to 
know how to sum up these logarithms. 

As a matter of fact the expression ( p. ID is not quite accurate theoretically. One 
should not forget that, in doing the loop integrations, in C„ we must discard the 
domain of virtual momenta below /i, by definition of Cn{fi)- Subtracting this domain 
from the perturbative loop integrals we introduce in C„ power corrections of the 
type (/i/mg)" by hand. In principle, one should recognize the existence of such 
corrections and try to learn how to deal with them. The fact that they are there 
was realized long ago (see e.g. V. Novikov et al, Ref. |^) and then largely ignored. 
If it is possible to choose fi sufficiently small these corrections may be insignificant 
numerically and can be omitted. This is what is actually done in practice. This is 
one of the elements of a simplification of the Wilsonean operator product expansion. 
The simplified version is called the practical version of OPE, see below. Certainly, 
at the modern stage of the theoretical development it is desirable to return to the 
issue to engineer a better procedure than just discarding these f^/mg terms in the 
coefficient functions. Attempts in this direction are under way [|6^] . 

Even if perturbation theory dominates in the coefficient functions they still con- 
tain also nonperturbative terms coming from short distances. Sometimes they are 
referred to as noncondensate nonperturbative terms. An example is provided by the 
so called direct instantons with the sizes of order mg . These contributions fall off 
as high powers of Aqcd/'tiq and are very poorly controllable theoretically. Since the 
fall off of the noncondensate nonperturbative corrections is extremely steep, basi- 
cally the only thing we need to know is a critical value of mg. For lower values of mg 
no reliable theoretical predictions are possible at present. For higher values of mg 
one can ignore the noncondensate nonperturbative contributions. There are good 
reasons to believe that the b quark, fortunately, lies above the critical point. Again, 
I must add that the noncondensate nonperturbative contributions are neglected in 
the practical version of OPE. 

(Do we see seeds of the nonperturbative contribution in Eq. (|5.1| )? Yes, we do. 
At any finite order the perturbative contribution is well-defined. At the same time, 
if the coefficients in the series ( |5.1| ) grow factorially with / - and this is actually the 
case - the tail of the series, / > l/a^, must be regularized which may bring in terms 
of order 

exp(-C/a.(mg))~ (^^y (5.2) 

where C is some positive constant and the exponent 7 need not be integer. In 
a sense, one may say that contributions to Cn of this type are vaguely related to 
diagrams with l/ag hard gluon loops.) 

Thus, two sources of nonperturbative corrections in the physical amplitudes are 
indicated. Those due to nonperturbative terms in the coefficient functions are sys- 
tematically ignored (and, perhaps, rightly so, as I tried to convince you) in these 



67 



lectures and in all works based on the practical version of OPE which constitute the 
overwhelming majority of all works devoted to the I/tuq expansions. The second 
source is operators of higher dimensions in the Lagrangian ( ^.11 ), the so called con- 
densate corrections. The latter were in the center of our attention; they generate 
the l/mq expansions discussed above. One new element which I would like to add 
here, is that the series of l/mg terms generated by higher- dimensional operators is 
also asymptotic and divergent in high orders Of course, we always calculate 



only one, at best two, first l/wig corrections, truncating the series. If, however, 
one would ask what the impact of the high-order tail of the power series is, the 
answer would be: this tail is reflected in exponentially small terms ~ exp(— mg). 
This type of contribution is certainly not seen in OPE truncated at any finite order. 
A transparent example is again provided by instantons. This time one has to fix 
the size of the instanton p by hand, po ~ A^^. Then their contribution to physical 
amplitudes is C(exp(— rngpo))- The relation between exp(— mgpo) piece and the 
high-order terms of the power series is conceptually akin to the connection between 
exp(— l/a^) terms and I ~ l/a^ orders in the perturbative expansion. 

Summarizing, Wilsonean OPE ( |5.1| ) leads to expansions in different parameters. 
Purely logarithmic terms (Inmg)"' are due to ordinary perturbation theory. Terms 
of the type (mg ) "'^ (In mg)^'^ reflect higher-dimension operators and direct instan- 
tons. In the former case the values of k are integer, the latter case may produce 
non-integer values of k. In the practical version of OPE we calculate the coefficient 
functions perturbatively. All non-perturbative terms come from condensates within 
this approximation. The condensate power series is truncated: only those operators 
whose dimension is smaller than some number are retained. 

The practical version of OPE was heavily used in connection with the QCD sum 
rule method. It was checked that in the majority of channels this is a valid 



approximation allowing one to calculate in the Euclidean domain down to p as low 
as 0.6 or 0.7 GeV. The validity of this approximation is an element of luck; it relies, 
among other things, on the fact that Aqqd is significantly smaller than 1 GeV, and 
as(lGeV)/7r is already a small parameter. 

I hasten to add that some exceptional channels where the practical version of 
OPE fails at much larger values of p were detected in the analysis of glueballs . 



It would be interesting to explore the issue in the context of the heavy quark theory. 
The existing theory gives no clues for establishing the domain of validity of the 
practical version of OPE from first principles, neither does it tell us about when the 
exponential terms, not visible by standard methods, become negligibly small. At 
this point we have to rely on indirect methods and phenomenological information. 

5.2 Untangling hard gluons from soft ones 



The coefficient functions C„ in Eq. (O) contain, generally speaking, an infinite 



number of perturbative terms, and non-perturbative contributions of different types. 
Practically we often calculate them to the first nontrivial order. For instance, in 
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Lecture 3 we treated the transition operator in the Born approximation; thus, all 
coefficients in OPE were found to order For a number of purposes (although not 
always, of course) such a calculation, ignoring the hard gluon exchanges altogether, 
is quite sufficient. Let me remind that by hard gluons I mean those with offshellness 
from n up to rriQ. Let us ask a question - can one find a theoretical parameter which 
would justify the approximation of no hard gluon exchanges? In other words, does 
a parameter exist that would allow one to switch the hard gluons on/off ? 
Each extra hard loop contains the running gauge coupling ^^(/i). 
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where b is the ffist coefficient of the Gell-Mann-Low function, 



(5.3) 



b = —N^ Uf. 

3 3^ 

If we could make b very large the running law of would be very steep effectively 
switching off all hard gluons. Indeed, once fi is bigger than, say, 2Aqcd and b oo 
the gauge coupling constant asifJ^) 0. The ffist idea which immediately comes to 
one's mind is to make b large by tending the number of colors Nc to infinity. Alas, 
this idea does not work. It is known from the early days of QCD that the expansion 
parameter in all planar diagrams is NcOg, not as itself p|. Thus, the diagram of 
Fig. 13 is of the same order in Nc as the Born graph of Fig. 1. So, we have to rely 
on numerical smallness of 1/6. For instance, in the theory with three light flavors 
and three colors 6 = 9, quite a large number. This is not the ffist time in physics 
we have to deal with numerical enhancements. It is true that it is always better to 
have an adjustable parameter, which could be sent to infinity at will, than to deal 
with just a large fixed number. It is quite unfortunate that we do not have such a 
parameter at our disposal in the real world QCD. If one still wants to have b as an 
adjustable parameter one could try a trick. Let us assume that, apart from quarks 
and gluons, our theory contain quark ghost fields. These ghost fields are perfectly 
the same as the quark fields, with a single exception - each ghost loop has an extra 
minus sign. The quark ghost fields may or may not have a mass term. Let us say 
that they do have a mass term nigh equal to Aqcd- Then they would automatically 
decouple in the soft contributions. The action of such a crazy theory has the form 

iS = iSqcB + ^9h{^ V> - ^gh)(lgh = 

q 

iSqcB - Ngh In det {i p - nigh} (5.4) 

where S'qcd is the action of quantum chromodynamics, see Eq. ( |l.lj ), and Ngh is 
the number of the quark ghosts, a free parameter assumed to be large. Notice the 
ghostly minus sign in front of the logarithm of the determinant. After some thinking 
one may conclude that, perhaps, this theory is not so crazy. Let us postulate that 
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the initial particles we consider belong to our world - B, D and so on - i.e. they 
do not carry these quark ghosts. Of course, if B decays the quark ghosts do appear 
in the final state, and the probability of their emission is negative. This does not 
mean, however, that the total amplitude is not unitary, as one could suspect from 
the fact that we introduced the fields with a wrong metric. Indeed, it is obvious 
that the only role of the quark ghosts is to switch off all hard gluons in the limit 
Ngh —>■ oo since in this limit 

b = ll-'^Nf + '^Ngh^oo 

and as{fi) 0, according to Eq. ( |5.3|) . In particular, the diagram of Fig. 13 where 
the gluon line is dressed with the bubble insertions vanishes. All soft contributions 
with /i < Aqcd remain intact, however, and the positivity of the forward scattering 
amplitudes is not violated. 

If there exists a stringy representation of QCD it should refer to the fake "QCD", 
Eq. (|5.2|) , rather than to the real one since in the string amplitudes there is no place 
for hard gluons. 

The idea of treating 6 as a numerically large parameter is not new in QCD. In 
the purely perturbative calculations it constitutes the basis of the so called BLM 
approach [^. Originally the BLM approach was engineered as a scale-setting pro- 



cedure intended as a substitute for full computations of O^a"^) corrections. Assume 
that 0{as) corrections in some amplitude are known exactly. In order typically 
one has to deal with a large number of graphs. The idea is to pick up only those 
which contain a "large parameter", presuming that the graphs without b are 
numerically suppressed. Typically there are very few graphs producing ba^. By do- 
ing so we can approximately determine the scale /i in the 0{as) term without labor 
and time-consuming calculation of a large number of all contributions. Later, 
it was suggested |]7T], to extend the prescription of the "6 graph dominance" to 



even higher orders, a more extremist and dangerous approach. In both cases the 
limit of large b is used to get some information about perturbation theory. I use this 
limit in order to switch off the perturbative hard gluons in the first place pushing 
the theory to the mode where only the soft gluons survive, hopefully providing a 
more transparent picture of the infrared dynamics determining the regularities of 
the hadronic world. 



5.3 Impact of hard gluons 

Having said all that let us return to the real world where b is fixed, not infinity, and 
examine several examples of corrections due to hard gluons. An instructive example 
to begin with is the calculation of the coefficient in front of the chromomagnetic 
operator Oq in the effective Lagrangian £heavy(/^), see Eq. ( |OD This coefficient 

^^In the limit 6 — > oo the coefficient given in this expression does indeed vanish, in fuU accord 
with the argument of the previous section. 
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takes into account virtual gluons with offshellness from n to mg. 

The hne of reasoning is as follows. Our starting point is n = mq. At this 
normalization point the Lagrangian we deal with is the QCD Lagrangian ( |1 . 1| ) with 
the coupling constant and heavy quark mass normalized at uiq. We then descend 
a little further, down to equal to a finite fraction of mg, say, mg/S. This is 
sufficient to make the Q quark nonrelativistic and make all nonrelativistic expansions 
work. Being interested only in logarithms of mq we ignore any nonlogarithmic 
corrections that may appear at this stage. The nonrelativistic expansion of the 
Lagrangian Q{i Jf) — raQ)Q implies that the operator Q{i/2)aGQ appears with the 
coefficient Cq = l/(2mQ). Further evolution down to ^ = several units xAqcd will 
change C; in particular, at one loop 




Co C{n) =Co{l+ 7-^ In — ^ + non-log terms (5.5) 



where 7 is a number. Our goal is to find 7 and then sum up all leading logarithms. 
This is not the end of the story, however, if one wants to represent the result in the 
form (|1.3|) , where the sum over the operators includes only the Lorentz invariant 
ones. The leading operator is ^heavyif^)- The coefficient C(/x) should be represented 

as 

C(/i) = Co + (Cifi) - Co) ; 

then Co is swallowed back in the definition of jC^^^-^y^fi) while the expression in the 
brackets represents Cg in Eq. ( |1.3| ). 

The relevant one-loop graphs are depicted on Fig. 14. At first sight the number 
of diagrams is rather large, and the computation might seem rather cumbersome. 
My task is to reduce it to a back-of-the-envelope calculation by using several smart 
observations and the background field technique. 

First of all, as it was already mentioned, we will be hunting only for the terms 
containing In mq / /x omitting all as terms without logarithms. The logarithms 
hamq/fi have a dual nature - they appear from the loop integrations where the 
integrands presents an infrared limit with respect to heavy quarks Q while presenting 
simultaneously the ultraviolet limit with respect to gluons. That is why they were 
called hybrid in Ref. |^|, the paper where these logarithms were discovered. In the 



language of HQET they are referred to as matching logarithms. 

Secondly, in this perturbative calculation we will naturally discard all l/mq 
corrections. 

The closed circle on the diagrams of Fig. 14 denotes the vertex {i/2)a^^Gij = 
—aB. Let us consider for definiteness only one term with i,j = 1,2, i.e. —a^Bz 
keeping in mind that other terms will give the same. 

It is absolutely obvious that the graph of Fig. 14e gives no contribution in our 
approximation. Indeed, the very existence of this graph is due to the nonlinear term 
[^1^2] in the definition of G12. However, neither Ai nor A2 interact with the heavy 
quark in the leading in l/mq approximation, as it is clear from Eq. ( |1.15| ), only Aq. 
(We work in the rest frame of the heavy quark Q.) 
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Next, let us analyze the diagrams c and d. To this end it is convenient to write 
the gluon Green function in the background field. For a detailed exposition of the 
technique the reader is referred to the review paper pS. For our purposes we need 



so little that it is quite in order to carry out all necessary derivations here. Let us 
split the four-potential in two parts - the external field (v4^)ext and the quantum 
part which will propagate in loops, 

= (v4^)ext + (5.6) 

As explained in Ref. the gauge conditions on (A^)ext and may be different, 
for instance, the Fock-Schwinger gauge with respect to the background field and 
the Feynman gauge with respect to the quantum field. Here we do not need to 
discuss the gauge condition on (A^)cxt- The quantum field will be treated in the 
Feynman gauge. The definition of the gluon propagator in the background field is 
standard: 

K = iTK{x),alm). (5.7) 
The Lagrangian of the quantum gluon field in the Feynman gauge has the form 

C = -^(/^r<)' + ga^iGlUatr'^ (5.8) 

plus cubic and higher order terms in plus the ghost terms - all irrelevant for the 
calculation at hand. Here 

D^J^X = d,al + (?/"''^(4)ext< . 



The second term in the Lagrangian ( |5.8| ) describes the interaction of the magnetic 
moment of the gluon quantum with the background field. If we switch off this 
magnetic terms for a short while we immediately observe that both graphs. Fig. 14c 
and (i, vanish. Indeed, the Lorentz structure of the first term in Eq. (|5.^ ) is such 
that the Green function generated by it is obviously proportional to g^v Hence the 
loops displayed on Figs, c and d can not be formed. Say, the diagram c requires 
converting the quantum leaving the vertex into the quantum coupled to the 
heavy quark. Let us now switch on the magnetic term and take into account the fact 
that the background field is chromomagnetic, not chromoelectric (I remind that we 
are interested in the vertex —aB.) This means that the graph c still vanishes since 
the conversion of Oj into can only take place in the chromoelectric background 
(Ggf*)- The diagram d is not vanishing, however, and is readily calculable. We start 
from the vertex 

make one insertion of the magnetic term in the Lagrangian (|5] 
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and after that one can take the free gluon propagators, which yields 



{^/2)2g'Qa''r'HVr''GU-^r^ J 1^ (5.9) 

where the factor 2 comes from two different ways of pairings. The integral over dk is 
evidently logarithmically divergent both at the upper and lower ends and should be 
cut off at TTiQ from above and at fi from below. This logarithmic divergence should 
be welcome since in this way we are going to get the desired hybrid logarithm. 
Equation ( |5.9| ) immediately leads to 

where Nc is the number of colors (A^^^ = 3). In other words the factor produced by 
the one-loop graph of Fig. 14d is 

7 iFig. i4d= -2Ar, . (5.10) 

The last step in our exercise is calculation of the diagrams of Fig. 14a and b. 
The Feynman integral for the diagram a is quite trivial, 

d^k ^ 1 



(5.11) 



where k is the virtual gluon momentum. Now, a minute reflection shows that in 
the Abelian theory (i.e. if the gluons were photons and the diagram of Fig. 14a 
was considered in QED) this contribution must be exactly canceled by that coming 
from diagrams b. This assertion can be traced back to the nonrenormalization of 
the Q'y^Q vertex in QED. (We should take into account the fact that the hybrid 



logarithms do not depend on the Lorentz structure of the vertex at all |£0[ and are 
the same for 7^ and a^jy.) This observation implies that in QCD the net effect of 
the two diagrams 146 reduces to replacing Eq. ( |5.11|) by 



Ncg' fd^k- ^ ^ i 1 

The ie prescription indicated explicitly defines the integration contour (Fig. 15). 
We first do the k^ integration using the residue theorem, then the remaining d^k 
integration and arrive at 

{^/2)Q{-a,Blt^)Q\n(^\ (5.13) 
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leading to the following "dressing" factor due to diagrams 146 and c: 



7 I Fig. Ua+b= Nc . . (5-14) 

The diagram 14/ must be discarded in the background field calculation - it 
merely renormalizes the gauge coupling constant included in the definition of Og- 
The overall one-loop dressing factor is obtained by adding up Eqs. ( |5.1(]| ) and 



7 = -iV, = -3. (5.15) 

Now the renormalization group allows us to sum up all leading log terms, in 
a standard manner; the summation leads to Eq. ( |1.4|) . The same result can be 
rephrased as follows: in the l/mg expanded effective Lagrangian jCheavy(/^) the over- 
all coefficient in front of Oq is 

This is nothing else than the reflection of the hybrid anomalous dimension of the 
operator Oq found in Ref. 0. 

It is curious to note that C^r, the second operator of dimension 5 (see Eq. (|3.6|) ), 
has vanishing anomalous dimension which can be proven with no calculations in no 
time. 

To see that this is indeed the case we merely repeat the argument preceding 
and following Eq. ( ^.5|) . Let us assume for a short while that the hybrid anomalous 
dimension of the operator (9^ is non- vanishing. Then after evolving to a low normal- 
ization point its coefficient gets renormalized, and there is no way one could absorb 
Ott back into a Lorentz invariant expression Q(J/ — mq) in the effective Lagrangian. 
Needless to say that £heavy(/^) (before the l/mg expansion) must be expressible in 
terms of the Lorentz invariant structures. 

A close line of reasoning leading to the same conclusion takes advantage of the 
expansion ( pl3|) , 

QQ - -^Q'-aGQ = QjoQ - t^Q^'Q + - (5.17) 

where the dots denote terms of the higher order in l/mq. The left-hand side is 
Lorentz scalar while the right-hand side is written as a sum of terms that are not 
Lorentz scalars individually. The matrix element of Q^oQ has the meaning of energy 
E (which at small velocities reduces to m + jf /2m) while that of the second term on 
the right-hand side has the meaning of —ff /2m. The first term is not renormalized 
by the gluon dressings, of course. If the coefficient of the second term was distorted 
by the anomalous dimension, the cancellation of the Lorentz noninvariant part would 
be ruined, and the right-hand side could not be equal to the left-hand side. 

Concluding this section let us discuss the impact of the hard gluons on the scaling 
law of, say, pseudoscalar coupling fp defined in Sect. 3.4. In this section it was 
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— 1/2 

shown that fp ~ mg modulo logarithmic corrections. Now we address the issue 
of the logarithmic corrections due to the hybrid anomalous dimension of the current 
Qliilbl- Let us add to the original QCD Lagrangian the term A£ = A^Q'-y^'-y^q 
where is an auxiliary c-number field and evolve A£ down to ^. The result of 
this evolution is the anomalous dimension 



the subscript here indicates the normalization point. The corresponding calculation 
is even simpler than that of the anomalous dimension of Oq and will not be discussed 
here. The interested reader is referred to Ref. or to review papers |Q Corre- 
spondingly the complete asymptotic scaling law of fp is fp ~ rriQ^^'^ {as{mQ))~'^^^ . 

5.4 fi dependence of the basic parameters of the heavy quark 
theory. Measuring A{fi) 



The Lagrangian ( p..3| ) summarizes the evolution from a high normalization point 



down to /i. Since all operators in this Lagrangian are normalized at fi it is perfectly 
natural that their matrix elements are also fi dependent. In particular, the matrix 
elements of Oc and denoted by /i^ and /i^ in Lecture 3 depend on /x. Actu- 
ally, Hq depends on fi rather strongly, through logarithms of /i - this is explicitly 
demonstrated by the fact that {OG)mQ = (q^s(/^)/q^s("^q))^'^^^(C^g)m- ^^is case 
hardly anybody would even think about tending ^ —* 0. As for the operator 0-,^, 
the situation here is trickier. As we saw, it has no diagonal anomalous dimension, 
still some /x dependence appears through mixing with QQ, see for details. 

Let us discuss now A, another basic parameter of the heavy quark theory. The 
issue of its fi dependence was at the epicenter of a heated debate recently. By itself A 
never appears in £heavy; moreover the quark mass mg appears in the l/mg expanded 
effective Lagrangian ( p. . 1 1|) only through l/mg corrections. Therefore, in the limit 
mg — > oo (which is often identified with HQET) it is quite tempting to say that 
A = Mhq — mg is a universal constant. For a few years it was taken for granted 
that such a constant exists. Within the framework of our approach based on the 
Wilsonean treatment of full QCD it is perfectly clear that this is not the case. The 



quark mass in Eq. (O) explicitly depends on /x resulting in a yU dependence of A. 



Since the issue is of importance let us rephrase this statement as follows. Since 
quarks are permanently confined the notion of the heavy quark mass becomes am- 
biguous. To eliminate this ambiguity one must explicitly specify the procedure of 
measuring "the heavy quark mass" . The definition through the effective Lagrangian 



^^Thc wording in these reviews is somewhat different. You will read about the matching log- 
arithms of HQET for the axial current. Technically this is perfectly the same as the anomalous 
dimension within our approach. 
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( p. .31 ) is consistent. Other definitions are certainly conceivable; any consistent proce- 
dure will necessarily involve a cut-off parameter /i , and then A(/i) = Mq — mqlfj,) 

The effective Lagrangian (|1.3|) is not something you directly measure, neither 
is rriQ. Defining A or rriQ is equivalent to saying how they are measured, how the 
parameters in the effective Lagrangian are related to measurable quantities. To this 
end one can use any suitable prediction of the heavy quark theory, in particular, 
Voloshin's sum rule ( |2.35| ). To avoid inessential technicalities I will discuss the issue 
in the framework of the toy model of Sect. 2.1. All results can be immediately 
extended to the real QCD. Equation ( p.35| ) gives a nice definition of A in terms of a 



measurable quantity, the average value of Eq'^^'^ — E where E is the energy of the 
quantum. The problem is that in Sect. 2.1 we discussed the question switching off 
all hard gluons, so that the above average value looked like a ^ independent number. 
To see where the /i dependence comes from we must include hard gluon corrections. 

If the gluon field is treated only as a soft medium the spectrum of the decay 
Hq ^ Xq + (p looks roughly as on Fig. 9. The shoulder to the left of the elastic 
peak arises due to production of the excited states. It is important that in this 
approximation the spectrum rapidly (exponentially) decreases outside the end point 
domain, so that the entire region of E from zero up to Eq^"^^— several units xAqcd 
remains unpopulated. The average 



/ ° dE^Ef^'-E) 
Jo 



1 dT 
f~odE 



is then a well-defined number independent of tuq or any cut-offs. 

The situation drastically changes once we include hard gluon emission. In cal- 
culating radiative gluon correction we can disregard, in the leading approximation, 
nonperturbative effects, like the difference between rriQ and Mhq or the motion of 
the initial quark inside Hq. Thus we deal with the decay of the free quark Q at rest 
into q + (l)+ gluon. The virtual gluon contribution merely renormalizes the constant 
h in the analysis presented above and is irrelevant. 

The effects from real gluon emission are most simply calculated in the Coulomb 
gauge, where only the graph shown in Fig. 13 contributes. A straightforward 
computation yields |2^, |3l|] to leading order in 



dV^'^ 8as E' 1 

dE ° 971 Eoml Eo-E' ^ ' 



^■^In the literature you can find assertions that an "absolute" heavy quark mass, or the so-called 
pole mass, can be defined and can be shown to be a universal number independent of any cut- 
offs. These assertions are false. The notion of the pole mass exists only to a given finite order of 
perturbation theory. No consistent definition of the pole mass can be given already at the level 
of the leading nonperturbative corrections ©(l/mg). The notion of the pole mass is absolutely 
foreign to the approach I present here, therefore, I do not want to go into details, see [Q. I will 
only say that it assumes it is possible to separate perturbative contributions from nonperturbative 
(?!) in contradiction with our approach which separates soft contributions from hard. 



76 



Hard gluon emission obviously contributes to the spectrum in the entire interval 
< E < Eq^^'^ creating a long "radiative" tail to the left of the end point domain, 
(note that in this calculation one can put Eq^^^ to Eq). In the first order calculation 
Us does not run, of course. Its scale dependence shows up only in the two-loop 
calculation; it is quite evident, however, that it is as^E — Eo) that enters. Therefore, 
strictly speaking, one cannot apply Eq. ( p.l8| ) too close to Eq. Even leaving aside 



the blowing up of as{E — Eq), there exists another reason not to use Eq. (|5.18 



m 



the vicinity of Eq: HE is close to Eq, the emitted gluon is soft; such gluons are to 
be treated as belonging to the soft gluon medium in order to avoid double counting. 

The separation between soft and hard gluons is achieved by explicitly introducing 
a normalization point /i. The value of /i should be large enough to justify a small 
value for as(/i). On the other hand we would like to choose /i as small as possible. 
We then draw a line: to the left of Eq — the gluon is considered to be hard, to the 
right soft. At E < Eq — fi the experimentally measured spectrum must follow the 
one-loop formula ( ^.18|) , see Fig. 16. 



Let us return now to Voloshin's sum rule, i.e. the first moment of Eq^^'^ — E, 
with radiative corrections included. A qualitative sketch of how dT/dE looks now 
is presented in Fig. 16. Because of the tail to the left of the end point domain we 
can not define A as the value of Eq'^^^ — E averaged over the entire range of the (p 
energy, < E < E^^^'^ . The integral would be proportional to agrriQ because of the 
domain of small E. Besides, this would contradict the physical meaning of what we 
want to define. By evolving the effective Lagrangian down to /i we include all gluons 
harder than /j, in uiq, thus excluding them from A. Thus, we must accept that 

A(/i)=/° ^^^-^i^Efy' -E)dE. (5.19) 
JEfy'-^^ v^TQdE^ ° ' ^ ' 

Since the explicit form of the tail to the left of the end point domain is known (for 
small as{fi)/TT the physical spectrum is supposed to tend to the perturbative result) 
the /i dependence of A becomes obvious, 

5A = 5^.^4^. (5.20) 
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Equation (|5.19|) provides us with one possible physical definition of A(yu) (among 



others) relating this quantity to an integral over a physically measurable spectral 
density. The pole-mass based definition, being applied to our example, would involve 
three steps: (i) Take the radiative perturbative tail to the left of the shoulder and 
extrapolate it all the way to the point E = Eq; (ii) subtract the result from the 
measured spectrum; (iii) integrate the difference over dE with the weight function 
{E — Eq^^^). The elastic peak drops out and the remaining integral is equal to 
Tq{v'^ /2)K. It is quite clear that this procedure cannot be carried out consistently 
- there exists no unambiguous way to extrapolate the perturbative tail too close to 
Eq^^'^ , the end point of the spectrum. Our procedure, with the normalization point 
introduced explicitly, is free from this ambiguity. 
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In practice, the fi dependence of A(^) may turn out to be rather weak. This is 
the case if the spectral density is such as shown in Fig. 16, where the contribution 
of the first excitations (lying within ~ Aqcd from E'q^^'^) is numerically much larger 
than the radiative tail representing high excitations. It is quite clear that if the 
physical spectral density resembles that of Fig. 16 and /i = several units xAqcd, 
the running A(yu) is rather insensitive to the particular choice of 

It remains to be added that a similar definition of A(/x) works in real QCD. Here 
it may be defined through an integral over the spectrum in the decay B — > X^lv 
measured in the domain where the recoil of the hadronic system is small, |^ <S mc, 
i.e. in the SV limit. 

5.5 Hard gluons and the line shape 

Actually we have already started considering this question in the previous section 
where the radiative correction to the spectrum in the decay Hq Xg + (f) was found 
in the SV limit. In this case the impact of the hard gluons is mild - they provide a 
long but squeezed tail outside the end point domain which could be evaluated in the 
leading (one-loop) approximation. The reason why there are no violent distortions 
of the spectrum is simple: the q quark produced is slow, and slow quarks do not 
like to emit hard gluons. If the final quark was fast it would produce gluons like 
crazy through bremsstrahlung, and the impact of such bremsstrahlung on the line 
shape would be much more drastic. As a matter of fact, if nig — > one can not 
limit oneself to any finite number of gluons - an infinite sequence of the so called 
Sudakov (or double-log) corrections must be summed over. 

By definition the Sudakov corrections are those in which each power of is 
accompanied by two powers of logarithm \nmQ/{mQ — 2E). When one approaches 
the end point domain the logarithm inevitably becomes large, and overcompensates 
the smallness of the gauge coupling constant ctg. So, the more gluons emitted the 
higher the probability. The phenomenon is classical in nature and has a transparent 
physical interpretation. Indeed, in the initial state Hq the color field in the light 
cloud corresponds to a static source. The final quark produced is very fast. The 
stationary state of the color field corresponding to a fast-moving color charge is 
strongly different from that of the stationary charge. Therefore, the excess of the 
color filed is just shaken off in the form of the multiple emission of gluons. If you 
forbid to emit a large number of gluons and insist that the final state is just "one 
quark" (this would correspond to the two-body decay kinematics and the delta- 
function-like narrow spectrum) then the probability of such an improbable event is 
terribly suppressed. This explains why for the massless final quarks the narrow peak 
in the end point domain obtained in Sect. 4.2 will be drastically distorted, and a 
well-developed tail to the left of the end point domain will appear. 

The theory of the Sudakov corrections constitutes a noticeable part of the per- 
turbative QCD, and here, of course, I have no possibility even to scratch the surface. 
I will give just a few hints referring the interested reader to the original papers and 
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textbooks |7|. 

The first order probability of emission of the massless gluon in the b ^ s'j decay 

is 

1 d^Tg 2a, 



TduM"^ 37rcj(l -cos^9) ^^'^^'^ 

where uj is the gluon momentum and d is its angle relative to the momentum of the 
q quark. In this expression it is assumed that uj -C mq. As a matter of fact it is 
perfectly legitimate to make this assumption since the double logarithm comes only 
from this domain of integration. The (j) energy in the presence of a gluon in the final 
state is given by 

^ ^ ml-2uj{mb-uj)(l-coiie) ^ tuq _ 

2mQ - 2u{l - cosO) 2 4a; ' ^ ' ^ 

k± ^ 00^ . 

One starts from computing the (first order) probability w{E) for the gluon to be 
emitted with such momentum that the quantum gets energy below given E. This 
probability is obtained by integrating the distribution ( |5.21 ) with the constraint 
that {mQ/2) — {kj_/4uj) is less than the given E, 

^E)= /d.d^^i/^4--f^-^)l=— In^^^- (5-23) 
^ ^ J Vdud^^ \4:U V 2 // Stt tuq - 2E ^ ' 

I integrated over d"^ first; the upper limit of integration is of order one, the lower 
limit is determined from the 9 function in Eq. ( |5.23| ). Then we can carry out the uj 



integration. The upper limit is mQ/2 while the lower limit is seen from the same 6 
function, 

u^<—-E. 

The function w{E) has the meaning of the probability of emission of a sufficiently 
hard gluon lowering the (f) energy below E. The all-order summation of double logs 
amounts then to merely exponentiating this probability [[74[ , 

S{E) = e-"'^^) . (5.24) 

The spectrum then takes the form 

idr_di 

TdE~ dE- ^^-^^^ 

We see that as E approaches the end point, E close to mQ/2, the spectrum gets 
suppressed, in full accord with our expectations, since the presence of the very hard 
(f) does not allow, purely kinematically, the gluon shower to develop and the color 
field to restructure itself. Notice that the Sudakov corrections merely redistribute 
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the probability, since the full integral over the spectrum remains unchanged. They 
pump events out from the end point domain to lower values of E. 

The double log approximation per se does not allow us to determine the scale 
of as in the Sudakov exponent S{E). For practical purposes the scale setting is 
of course very important since exp{—w{E)) is a steep function. The question goes 
far beyond the scope of this lecture. Some partial answers can be found in Refs. 
Il75| , [76[| , see also suffice it to mention here that in Eq. ( |5.23|) turns out 



to be [|76[ as{\j{rnQ — 2E)mQ). Another point deserving stressing is that with the 
classical Sudakov formula one can not travel over the energy axis too close to the 
end point E = mQ/2 (even after the scale setting). Indeed, if E > (mQ/2) — fi the 
gluons emitted become too soft; such gluons constitute the soft gluon medium and 
have nothing to do with the perturbative calculation; they have to be referred to 
the primordial distribution function. Equation ( |5.23|) is applicable provided that 



A < mg - 2E < mg . 

If we come closer to the end point domain the classical Sudakov factor must be 
modified by cutting off and discarding the contribution of the soft gluons. This idea 
gained recognition only recently; it is obviously premature to further immerse into 
this topic for the time being. 

If the effect of the hard (perturbative) gluon emission is known the full physi- 
cal spectrum is obtained by convoluting the perturbative one with the primordial 
distribution function, 

-dE- = ^^^^ J "^y^^y^ dE • ^'-2') 

Integration over y runs from — oo to 1 (more exactly, the lower limit of integration 
is yo = — mg/A but this difference can be ignored). One should keep in mind that 
dV^Q^^/dE is nonvanishing only in the interval (0,mg/2). The convolution formula 
above is legitimate only as long as one does not apply it to the very low energy part, 
~ A; further details are presented in Ref. . 
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6 Figure Captions 



Fig.l. The forward scattering amplitude Q ^ Q. The dashed hne denotes the 
quantum. The sohd hne connecting two vertices is the q quark Green function 
in the background gluon field. The thick solid lines describe the Q quarks in the 
background field. 

Fig. 2. The transition operator relevant to the total semileptonic width of the 
heavy mesons. The dashed lines denote the leptons, / and u. Other notations are 
the same as on Fig. 1. 

Fig. 3. The two-point function ( p.63| ). The wavy line is the external scalar or 
pseudoscalar current; the quark propagator is in the background field. 

Fig. 4. The two-point function of two scalar currents in the scalar QCD. 



Fig. 5. The four-point function appearing in Eq. (|2.64|) . The dashed line denotes 
the current Q^Ga/sQ- 

Fig. 6. The three-point function relevant to the proof of the Isgur-Wise formula. 

Fig. 7. The "photon" spectrum in the decay Hq Xqcj). The final quark is 
assumed to be massless. The thick line represents the delta-function spectrum of the 
free quark approximation. The solid line is a sketch of the actual hadronic spectrum 
in the end point domain (possible radiation of hard gluons is neglected). 

Fig. 8. Evolution of the spectrum of Fig. 7 as the mass of the final quark 
increases (the schematic plot refers to iriq ~ mqll. The effects of the hard gluon 
bremstrahhlung are not included. 

Fig. 9. The photon spectrum in the SV limit, t;^ ^ 1. The dashed line shows 
the would be spectrum of the free quark decay. 



Fig. 10. Kinematically allowed domain in the transition B — > Xylv. The thick 
line indicates the populated phase space in the free quark decay. The shaded area 
of width ~ A is the end point domain populated due to the primordial motion of 
the 6 quark inside B. The shaded square in the lower right corner is the exclusive 
resonance domain where the inclusive approach developed here is inapplicable. 



Fig. 11. More or less realistic light-cone primordial distribution function versus 
X (borrowed from Ref. ||63|| ). 



Fig. 12. Kinematically allowed domain in the decay B —>■ XJ,v. Two large 
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circles show the domains where description should be based on the light cone and 
temporal distribution functions, respectively. 

Fig. 13. Correction of the first order in ctg to the transition operator of Fig. 1. 
Shown is the only graph contributing to the imaginary part in the Coulomb gauge. 
The gluon is denoted by the curly line. 

Fig. 14. One-loop diagrams determining the coefficient cq- The wavy hne 
denotes the gluon quanta, dashed line background glion field. 

Fig. 15. Integration contour in the plane. 

Fig. 16. A sketch of the photon spectrum in the SV hmit with the hard gluon 
radiation included. 
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